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I.  INTRODUCTION 


This  is  tlie  final  report  for  Contract  N00014-86-C-2055.  It  covers  research  performed 
by  Berkeley  Scholars,  Inc.  from  6  December  1985  to  5  December  1987;  ‘  The  research 
involved  the  theoretical  and  numerical  analysis  of  the  physics  of  free  electron  lasers  using 
relativistic  particle  bemns. 

The  major  emphasis  of  the  research  was  to  obtain  design  criteria  for  the  development 
of  a  two-stage  FEL  oscillator  operating  in  the  trapped  particle  mode,  and  for  a  UV-FEL 
operating  in  an  oscillator  configuration  with  an  intense  CO2  laser  beam  as  a  pumping 
source  instead  of  the  usual  wiggler  field.  In  order  to  carry  out  this  program  it  is  necessary 
to  provide  an  appropriate  model  for  the  radiation  physics  of  both  types  of  FEL’s,  to 
examine  the  relevant  nonlinear  wave-particle  dynamics,  and  to  model  the  characteristics, 
development  and  evolution  of  the  associated  electron  beam. 

Task  1,  to  develop  a  general  model  for  the  radiation  physics  occurring  in  a  FEL  with 
specific  application  to  a  two-stage  FEL  oscillator  operating  in  the  trapped-particle  mode, 
and  to  a  UV-FEL  operating  in  an  an  oscillator  configuration  that  uses  an  intense  CO2 
laser  beam  as  a  pumping  source,  is  a  necessary  prerequisite  for  Task  2  and  Task  3.  These 
tasks  are  model  applications,  and  the  research  addressing  these  tasks  has  given  rise  to  new 
results.  This  final  report  documents  the  research,  and  the  six  appendices  contain  the  details 
of  the  work.  Articles  prepared  from  the  appendices  have  been  submitted  to  professional 
journals  or  technical  conferences.  These  articles  have  been  prepared  in  collaboration  with 
scientists  from  the  Plasma  Theory  Branch  at  the  Naval  Research  Laboratory. 

The  effects  of  radiation  damping  on  beam  quality  in  the  inverse  free  electron  laser 
accelerator  are  discussed  in  Appendix  I,  also  published  in  the  journal  Particle  Accelerators. 
A  beam  envelope  equation  is  derived  and  solved  for  an  arbitrarily  tapered  wiggler  field. 
The  expression  for  the  evolution  of  the  normalized  transverse  beam  emittance  is  derived 
and  found  to  decrease  exponentially  with  distance  due  to  radiation  damping  until  it  is 
limited  by  quantum  excitation.  Our  results  show  that  substantial  improvements  in  beam 
quality  can  be  realized  for  acceleration  distances  comparable  to  the  radiation  damping 
e-folding  length. 

The  subsequent  appendices  give  an  account  of  the  subsequent  part  of  our  research, 
viz.,  analy.ses  directly  addressing  the  electron  beam  quality  issues  that  form  task  three  of 
the  statement  of  work.  Appendix  II,  also  published  in  Physical  Review  Letters,  concerns 
the  nonlinear  particle-wave  proce.sses  that  can  give  rise  to  radiation  focusing  and  guiding 
in  a  free  electron  laser.  If  the  centroid  of  the  electron  beam  is  trans\-ersely  displaced  the 
radiation  can  be  guided  by  the  electron  beam.  A  spatial  modidation  on  the  electron  beam 
envelope  can  also  modify  the  radiation  field.  These  and  other  phenomena  are  studied  using 
a  novel  source  (k'pendent  modal  representation  of  the  fully  three  dimensional  radiation  field, 
the  Source  Depc'iident  Expansion  (SDE)  method.  Among  the  merits  of  this  approach  is 
that  few  modes  an'  needc'd  for  an  accurate  description  of  the  radiation. 
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Appendix  III,  also  published  in  the  Physical  Review  A,  is  a  further  study  of  the 
focusing  and  guidance  of  the  radiation  by  the  electron  beam  using  the  SDE  method. 
Fast  and  accurate  numerical  solutions  of  the  fully  three  dimensional  FEL  problem  can 
be  obtained  over  distances  of  many  Raleigh  lengths.  The  effects  of  finite  einittance  and 
wiggling  of  the  electron  beam  can  also  be  studied.  Appendix  IV,  presented  at  the  1986 
International  FEL  conference  in  Scotland,  treats  similar  material.  Appendix  V,  presented 
at  the  1986  Conference  of  Particle  Accelerators,  expands  on  this  work  by  examining  the 
perfectly  guided  radiation  beams  in  the  Compton  exponential  gain  regime  of  an  FEL. 

Appendix  VI,  also  published  in  the  Physical  Review  A,  uses  the  SDE  method  for 
a  related  problem,  viz.,  the  optical  gain,  phase  shift,  and  spatial  profile  of  the  coupled 
electron  and  radiation  beams  in  an  FEL.  This  research  was  carried  out  in  collaboration 
with  B.  Hafizi  from  SAIC  and  with  NRL  personnel.  The  gain,  phase  shift,  wavefront 
curvature  and  radius  of  the  radiation  envelope  in  a  free  electron  laser  amplifier  are  obtained 
in  the  small  signal  regime.  The  electron  beam  is  assumed  to  have  a  Gaussian  density 
distribution  in  the  transverse  direction.  Numerical  calculations  indicate  that  the  radius  and 
curvature  of  the  radiation  beam  entering  a  wiggler  asymptote  to  unique,  spatially  constant 
values  after  a  finite  transition  region.  However,  in  the  asymptotic  regime  the  wavefronts 
diverge.  Analytical  expressions  for  the  gain,  phase  shift,  curvature  and  spot  size  axe 
derived.  It  is  shown  analytically  that  small  perturbations  of  the  optical  waist  and  curvature 
about  the  matched  value  are  spatially  damped  out,  indicating  the  stability  of  the  matched 
envelope.  When  the  electron  beam  is  modulated  in  speice,  the  optical  spot  size  oscillates 
with  an  almost  identical  wavelength  but  delayed  in  phase.  In  the  case  of  small  eimplitude, 
long-wavelength  betatron  oscillation  of  the  electron  beam  envelope,  generation  of  optical 
sidebands  in  wave  number  space  is  examined.  The  elTect  of  the  dispersion  characteristics 
of  the  primary  wave  is  found  to  be  negligible  for  typical  experimental  parameters. 
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II.  TECHNICAL  DISCUSSION 


Electron  beam  quality  as  measured  by  the  transverse  emittance  is  usually  determined 
by  the  gun  and  propagation  configurations  in  accelerators  (see  Ting  and  Sprangle,  1987 
for  this  discussion).  Under  idealized  conditions,  the  transverse  normalized  beam  emittance 
remains  a  constant  of  motion  as  the  beam  propagates  through  the  accelerator.  Therefore, 
to  improve  the  quality  of  the  beam,  it  is  necessary  to  decrease  the  beam  emittance  at  the 
injection  point.  However,  since  the  normalized  beam  emittance  is  essentially  the  trans¬ 
verse  area  in  phase  space  for  the  collection  of  beam  particles,  one  can  in  principle  reduce 
the  emittance  if  a  dissipative  mechanism  is  introduced.  A  natural  candidate  for  such  a 
dissipation  mechanism  is  the  induced  synchrotron  radiation  damping  due  to  the  transverse 
motion  of  the  particles  in  an  external  periodic  transverse  magnetic  field.  It  is  this  mecha¬ 
nism  that  will  be  focused  on  in  Appendix  I  when  the  external  magnetic  field  is  chosen  to 
be  a  helical  wiggler  field.  Since  this  radiation  damping  effect  is  small  at  low  energies,  it  is 
in  the  context  of  the  recently  proposed  high  energy  IFEL  accelerators  (references  1-11  in 
Appendix  I)  that  will  be  emphasized  and  concentrated  on  in  Appendix  I. 

First,  the  electron  orbits  in  an  IFEL  accelerator  must  be  obtained.  A  fully  relativistic 
formulation  of  the  equations  of  motion  which  include  radiation  damping  force  is  consid¬ 
ered.  The  damping  coefficients  are  obtained  from  the  transverse  dynamics  of  the  particles 
while  the  axial  dynamics  describe  the  acceleration  of  the  particles.  In  the  second  section, 
a  relativistic  envelope  equation  for  the  average  radius  of  the  electron  beam  is  derived, 
assuming  continuous  emission  of  the  synchrotron  radiation.  It  is  apparent  from  this  en¬ 
velope  equation  that  the  normalized  transverse  emittance  decays  exponentially  at  a  rate 
given  by  the  radiation  damping  coefficient.  The  envelope  equation  is  solved  using  a  WKB 
method  and  the  spatial  evolution  of  the  beam  radius  is  obtained.  Quantum  excitation  sets 
a  minimum  value  on  the  normalized  transverse  emittance  in  an  IFEL  accelerator  and  it  is 
derived  in  the  fourth  section  of  Appendix  I.  Strong  focusing  is  found  to  be  necessary  to 
reduce  such  minimum  to  an  acceptable  value.  An  example  is  given  in  the  last  section  of 
Appendix  I  for  a  set  of  proposed  IFEL  accelerator  parameters  (reference  2  in  Appendix  I). 
It  is  found  that  radiation  damping  does  reduce  the  emittance  of  the  accelerated  electron 
beam  while  resulting  in  an  insignificant  loss  in  particle  energy. 

In  analyzing  radiation  focusing  and  steering  in  the  FEL  by  using  a  source  depen¬ 
dent  expansion  technique,  it  is  found  that  in  the  one-dimensional  analysis  of  the  FEL 
the  radiation  field,  wiggler  field  and  electron  beam  resonantly  couple  so  as  to  modify  the 
longitudinal  wave  number  of  the  radiation  field  (references  1-3  in  Appendix  III;  also  see 
Sprangle,  Ting  and  Tang,  1987a  for  this  discussion).  This  can  lead  to  focusing  of  the 
radiation  beam,  a  phenomena  which  has  been  shown  to  play  a  central  role  in  the  practical 
utilization  of  the  FEL  (reference  9  in  Appendix  III)  since  in  many  proposed  experiments 
the  short  wavelength  radiation  beam  will  not  be  confined  or  guidc'd  by  a  waveguide  struc¬ 
ture.  Furthermore,  the  interaction  length  (wiggler  length)  is  usually  long  compared  to  the 
Rayleigh  length  associated  with  the  radiation  beam.  Therefore,  focusing  of  the  radiation 
beam  via  the  resonant  interaction  with  the  electron  beam  is  riect'ssary  in  order  to  ov('rcom(' 
the  natural  tendency  of  the  radiation  beam  to  diffract.  If  diffraction  of  the  radiat  ion  field 
were  not  fully  or  i)artially  offset  i)y  the  focusing  effect,  the  FEL  would  suffer  from  rechicecl 
gain  and  efficiency. 
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Two  primary  objectives  are  associated  with  the  use  of  a  source  dependent  expansion 
technique  in  the  analyses  presented  in  Appendix  III  (Sprangle,  Ting  and  Tang,  1987a)  and 
Appendix  IV  (Sprangle,  Ting  and  Tang,  lOSTb),  namely  to  present  a  general  method  of 
formulating  and  solving  problems  involving  radiation  focusing  and  guiding  for  in<'chanisms 
in  which  the  refractive  imlex  is  known  ;uid  to  apply  this  approach  to  the  focusing  and  steer¬ 
ing  of  radiation  in  FELs  with  mbitrary  gain.  The  method  is  :v  general,  self-consistent,  fully 
nonlinear,  modal  representation  with  application  to  the  phenomena  of  radiation  focusing 
and  guiding  in  FELs.  The  .source  function  (driving  current)  is  incorporated  self-consistent ly 
into  the  functional  dependence  of  the  radiation  waist,  the  radiation  wave  front  cur\'atvne 
and  the  radiation  complex  amplitude.  The  fundamental  mode  remains  dominant  through¬ 
out  the  evolution  of  the  radiation  field  due  to  the  source  dependent  nature  of  this  mo< 
expansion.  This  source  dependent  expansion  (SDE)  scheme  appears  to  have  a  num  ler 
of  advantages  over  the  conventional  vacuum  representation  (reference  7  in  Appendix  IV). 
Among  the  advantages  are  that  relatively  few  modes  are  needed,  compared  to  the  vacuum 
expansion  approach,  to  accurately  describe  the  radiation  beam.  Because  far  fewer  modes 
are  needed,  fast  numerical  solutions  of  the  fully  three-dimensional  wave  equation  can  be 
obtained  over  long  propagation  distances.  It  therefore  appears  feasible,  using  the  SDE  ap¬ 
proach,  to  incorporate  simultaneously  into  the  model  for  the  driving  current  density,  the 
effects  of  electron  beam  emittance,  energy  spread,  wiggler  gradients,  sideband  frequencies, 
etc.  Furthermore,  since  the  lowest  order  terms  in  this  expansion  are  a  good  approximation 
to  the  radiation  field  even  for  propagation  distances  long  compared  to  a  Rayleigh  length, 
valuable  insight  concerning  focusing  and  guiding  can  be  obtained. 
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I.  single  Particle  Dynamics 


We  shall  consider  the  motion  of  an  electron  under  the  Influence  of  a 
helical  wlggler  field  and  a  circularly  polarized  electromagnetic  wave  with  the 
inclusion  of  the  radiation  reaction  force.  The  fully  relativistic  equation  of 

1  ii 

motion  is 


^2.  V  X  B 

- |e|(E  ♦  ^=—)  *  , 


dt 


(1) 


where 


^R‘  2  2 
m  c 
o 


1  ‘*IeI 

r  ( V  )  ( - ) 

Y  dt  ^ 


2  dp  o  H  dp 


Is  the  radiation  damping  force,  ■  2jej^/5nQC^,  and  Y^  -  1  +  l^l^/m^c^. 

A  A 

The  radiation  field  is  given  by  Its  vector  potential  A  -A  (cos^e  -  sin^  ), 

“L  LX  y 

where  4  ■  kz-ut.  We  shall  assume  z-dependence  for  both  the  magnitude  and 

period  of  the  wlggler  field.  The  vector  potential  A  tor  the  helical  wlggler 

field  is  given  by  A  »A  [cosee  ♦sin6e  ]  where  A  «=  A  (z)  and  e  =  I  k.(z')dz''. 

°  — ww  x  y  ww  J  w 

^  o 

The  requirement  that  the  wlggler  field  satisfies  both  V  •  -  0  and 

V  X  B^  -  0  Introduces  transverse  variation  as  well  as  a  nonzero  z-oomponent  of 
the  magnetic  field. 

Since  we  shall  be  primarily  Interested  in  laser  driven  accelerators,  the 

normalized  wlggler  field  strength  a^  -  felA^/m^c^  is  assumed  to  be  much 

greater  than  the  corresponding  quantity  a^^  •  \e\k^^/m^o^  ter  the  radiation, 

i.e.,  a  »  a  .  It  can  then  be  shown  that  the  major  contribution  to  the 
W  L 

radiation  damping  Is  from  the  wlggler  field. 

We  shall  first  look  at  the  radiation  damping  term  In  Eq.  (1).  By 
neglecting  the  transverse  dependence  of  the  wlggler  field  for  a  beam  that  Is 
confined  sufficiently  close  to  the  axis,  we  have  the  Immediate  consequence 
that  the  canonical  momenta  in  the  x  and  y  directions  are  constants  of  motion 
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and  cay  be  chosen  to  be  zero.  The  mechanical  momenta  are  then  given  by 
I0 I  ^  I0  I  ^ 

Pjj  -  At*®x’  ^y  "  o'  -T*®y’  *  -L*  zeroth 

order  approximation,  the  total  relativistic  energy  is  conserved  which  leads 

to  Y  -  "  and  p  -  0.  Therefore,  the  only  significant  term  remaining  in  the 
z 

radiation  reaction  force  is 


F 


R 


.2 

d  £  ^  £ 

2  ~  22 
dt  me 


d£  2 


Neglecting  terms  of  order  «  1»  the  components  of  the  radiation  reaction 

R  R  H 

force  are  F  =  -  vicp  ,  F  -  -  vicp  ,  F  »  -  V|Cp  ,  where 
X  ‘^x’  y  I  ^y’  z  I  *'z' 


1 

"I 


^1.2  2 


(2a) 


(2b) 


are  respectively  the  spatial  decay  coefficients  due  to  radiation  damping  in 

2 

the  transverse  and  axial  directions.  Note  that  vi  -  vi  for  a  >>  1  which  is 

1  I  « 

the  case  in  the  IFEL  accelerator. 

The  most  significant  feature  of  the  transverse  motions  of  the  electrons 
is  the  betatron  oscillation  caused  by  either  the  Inhomogeneity  of  the  wlggler 
field  in  the  transverse  plane  or  other  focusing  mechanisms.  It  can  be  shown 
that,  for  small  oscillations  about  the  axis  of  the  wiggler  field,  the 
transverse  equations  of  motion  are. 


'  T  *  Y  ’ 


T  *  '’i’  ^  • 


(3a) 

(3b) 


where  d/dt  “  v  3/3z,  v  *  c,  *  ■  3/3z  have  been  used,  and  Kn  is  the  wave 
z  z  o 
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number  of  the  longitudinal  betatron  oscillation.  For  betatron  oscillations 

that  are  originated  from  the  v  x  B  force  due  to  the  nonzero  magnetic  field  in 

the  z-direction  of  the  realizable  wiggler  field, -  a  k  /{/zy)  , 

B  w  w 

The  axial  motion  of  the  electron  is  governed  by 


ym 


dv 

_ 7 

dz 


mv 


z  dz 


(v  X  B) 

z 


(H) 


where 


dz 


■|®ll  •  I  *  '>|P|  ^ 

3  2  3 

m  O'*  a  c^y 

o  o 


It  is  straightforward  to  show  that  the  axial  electron  acceleration  is 


2Y  3z 


3a  2a  a,  k  c  2k  - 

w  w  L  w  .  w  -!•  L  w 

- = —  sin  -  -  '^l^z  *  - 5 —  ’ 


(5) 


z 

where  +  J  Ck+k^Cz')  -  «/v^(z' )  jdz'  +  ip  is  the  phase  between  the 

o 

electrons  and  the  ponderomotive  wave  generated  by  the  beating  between  the 
radiation  and  wiggler  fields,  and  is  the  initial  phase  at  the  entrance  of 
the  interaction  region.  Equation  (5)  can  be  transformed  into  the  following 
pendulum  equation 


aS  k  ,  2v,k  3via  a  > 

— X  «  - -  _ —  ♦  - - -  sin  ip - • —  +  - ^ — 


dz 


dz 


2Y 


2  3z 


cos  <p 


(6) 


CY 


The  rate  of  change  of  relativistic  energy  may  be  obtained  from  Eq.  (H)  and  is 


a  a  k  i, 

_  -  sin  ♦  -  v.Y  ♦  -  2)  cos  i,  -  -^(a^ 

i  W 


dY 


Jr  2  2. 

- *-(a  ♦  a  )  . 

Lj 


(7) 
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Equations  (3).  (6)  and  (7)  will  be  the  basic  equations  we  shall  use  In 

studying  the  effects  on  beam  quality  due  to  radiation  damping.  The  terms 

containing  cos  li^  In  Eqs.  (5),  (6)  and  (7)  as  well  as  the  last  term  In  Eq.  (7) 

2  2  2 

may  be  neglected  when  the  conditions  a  »  a,  ,  a  >>  1,  k  >>  k  ,  and 

W  L  w  .  w 

2 

Y  >>  1  are  satisfied.  These  conditions  are  easily  achieved  in  high  energy 


IFEL  accelerators. 


II.  Derivation  of  Envelope  Equation  with  Radiation  Dampln 


The  single  particle  equations  of  motion  that  we  have  developed  in  the 
last  section  will  enable  us  to  study  the  macroscopic  behavior  of  the  beam. 
This  is  accomplished  by  considering  the  evolution  of  various  averaged 
quantities  associated  with  the  single  particle  variables.  *  We  begin  by 
multiplying  Eq.  (3a)  by  x*  and  x,  and  Eq.  (Sb)  by  y*  and  y,  where  ' 
denotes  3/dz.  Combining  the  resulting  equations  yields  the  following  set  of 
equations 


1  i.  fi2  ^ 

2  dz  2  dz 


2^2 

dz 


2  2  2 


-  Ji^p2 


2  dz 


dl  , 
dl  “  ' 


where  r^  -  x^  ♦  y^,  8^  -  x'^  ♦  y*^,  p  -  Y’/Y  ♦  vj^,  and  I  ■  (x'y  -  y’x)  is  the 

2 

normalized  angular  mom«itum.  We  eliminated  by  substituting  Eq.  (8b)  into 
Eq.  (Ba).  By  taking  transverse  ensemble  averages  over  beam  particles  in  Eq. 
(8),  and  denoting  the  ensemble  average  of  r^  by  a^  >  <r2>,  we  obtain  an 
equation  which  governs  the  evolution  of  the  root-mean-square  radius  of  the 
electron  beam, 

2d  2  d^2-„22d  ,u  d  2.  ^  1  d^a^  ^  d  ,„2  2. 

d^  ®  ♦  2pK  a  ♦  _(H  ^  )  4  -  _  4  — (K  a  ) 

dz  dz"' 

(9) 


It  is  easy  to  show  that  the  integration  factor  for  Eq.  (9)  is  g  a  where 

2  2  f* 

g  -  Y  exp(2j  vidz’)  .  Equation  (9)  can  now  be  put  into  the  form 
d/dz[  g^(a^a**  ♦  ya^a’  ♦  a  Kf)]  -  0  ,  and  can  be  integrated  to  give 

D 
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2  3  3  ^22  9 

g  [a-^a"  *  pa-^a’  ♦  a  Kg]  -  H  ,  where  Is  a  constant  of  znotlon  associated 

with  the  beam.  It  can  be  shown  that,  using  the  following  representation  for 

the  particles'  normalized  transverse  velocities,'^ 


a'  *  Lr  * 

■  r  'r  *  -2  *e  '  “i  • 


Where  la  the  normalized  transverse  velocity  spread,  and  L 
Eq.  (8c),  the  constant  is  given  by 


■  <l>  from 


-  Y^(0)L^(0)  +  Y^a^<|6e|l^>exp(2j  v^dz' ) 


where  Y(0)  *  Y(z-O)  and  L(0)  -  L(z«0) .  We  may  therefore  define  the  squared 

2  22  2 

normalized  beam  emittance  as  ■  Y  a  <|6ej^|  >  and  arrive  at  the 


following  envelope  equation 


d  a  .  1  dY  .  X  da  „2_ 

;;7  • '  i  31  *  I’  *  S" 


The  spatial  dependence  of  the  normalized  emlttance  and  average  angular 
momentum  are  given  respectively  by 


z 

c^(z)  -  c^(0)exp(-/  Vjdz*) 

n  -L. 


(11a) 


z 

L(z)  -  (Y(0)/Y)  L(0)exp(-|  vidz’) 

n  i 


(11b) 


where  e^(0)  *  c^(z=0).  Equation  (10)  together  with  Eq.  (11a, b)  constitute  the 
beam  envelope  equation  with  radiation  damping  terms  Included. 

One  can  see  that  when  vi  ■  0,  Eq.  (11a)  shows  that  e  remains  constant 

i  n 

and  Eq.  (10)  reduces  to  the  usual  relativistic  beam  envelope  equation  where 


8 


Is  the  familiar  normalized  beam  emittance.^^' Therefore,  in  the  presence 

of  radiation  damping,  the  root-mean-square  beam  radius  is  still  described  by 

an  envelope  equation  but  the  normalized  beam  emittance  is  no  longer  constant 

but  decays  exponentially  according  to  Eq.  (11a).  However,  the  decay  of  the 

normalized  beam  emittance  will  eventually  be  limited  by  quantum  excitation  due 

to  the  discrete  nature  of  the  synchrotron  radiation.  It  is  shown  In  a  later 

section  that  when  an  equilibrium  is  reached  between  these  two  competing 

processes,  the  minimun  normalized  emittance  achievable  through  radiation 

damping  in  the  IFEL  accelerator  is  given  by  (e^)  ,  «3fia^k  /(/Sn  cK,»)  . 

n  min  w  w  o  B 

In  the  presence  of  radiation  damping,  the  average  angular  momentum  also 
decays  exponentially  as  given  by  Eq.  (11b).  However,  one  may  choose  L(0)  »  0 
for  beam  generation  schemes  that  do  not  impart  an  average  angular  momentum  to 
the  electron  beam,  i.e.,  zero  magnetic  field  at  the  cathode.  We  shall  assume 
that  this  is  the  case  in  our  study  of  beam  quality.  He  shall  also  not 


,  and  will  denote  both  by  v. 
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III.  Evolution  of  Beam  Radius 


The  equation  for  the  root-mean-square  radius  a  In  Eq.  (10)  is  nonlinear. 
It  is  found,  however,  that  the  mean  square  radius  a^  satisfies  Eq.  (9),  which 
Is  a  linear  differential  equation.  For  beam  focusing  provided  by  the  wiggler, 
Eq.  (9)  may  be  solved  exactly  for  untapered  wiggler  fields  when  Y'  -  0. 

If  Y'  •>  0  or  when  the  tapering  Is  known,  it  can  be  solved  using  a  WKB  method 
if  we  assume  the  coefficients  are  slowly  varying.  Equation  (9)  can  be 
simplified  in  certain  limits  of  accelerator  designs  to  facilitate  analytical 
study.  It  can  be  shown  that,  Y’/Y  «  Kg  and  v  i  Kg,  which  allow  us  to  arrive 
at  the  following  approximate  equation 

S'"  +  3yS"  ♦  »»k|s'  ♦  C»»pKg  ♦2(k|)’]S  -  0  ,  (12) 


where  S  =  a^. 

In  order  to  obtain  net  acceleration  of  the  electrons  trapped  In  the 
ponderomotive  potential,  the  wiggler  field  must  be  spatially  tapered.  In  such 
a  case,  the  envelope  equation,  Eq.  (12),  is  a  linear  differential  equation 
with  spatially  dependent  coefficients.  We  solved  it  by  using  the  WKB-method 
which  assumes  these  coeffients  to  be  slowly-varying  functions  of  longitudinal 
distance.  By  assuming  both  K^/Kg  and  y  «  Kg,  the  general  solution  to 
Eq.  (12)  Is  found  to  be 

-M 

S  -  e  r^T-r  [  A  ♦  B  cos  2E  ♦  C  sin  2E  ]  , 

Kb(z) 

z  z 

where  M  -  J  y(z'' )dz' ,and  I  “  J  Kg(z')dz'.  The  coefficients  A,  B,  C  can  be 
o  o 

found  by  using  the  Initial  conditions  for  a  matched  beam,  a(z«0)  ■  a^, 

f  f  f 

a  (z-0)  -  0,  a  (z-0)  -  0  .  The  matched  beam  radius  a^  Is  related  to  the 

10 


initial  transverse  emittance  a 

o 


-  €^(0)/(Ko(0)Y^(0))  .  Using  the  initial 
n  D 

conditions,  we  arrive  at  the  following  expression  for  the  root-mean-square 
beam  radius. 


-M/2  t''B*®*i1/2i  , 

•  -  “o  *  ‘it^J  I  ’ 


-(0)  ♦  „i,/2 

2K.(0) 

D 


(13) 


Equation  (13)  shows  that  the  beam  radius  does  not  remain  constant  even  when 
the  beam  Is  matched  at  Injection.  In  addition  to  the  exponential  decay  from 
the  radiation  damping,  the  beam  envelope  developes  an  Induced  betatron 
oscillation.  However,  the  normalized  emittance  Is  Just  an  exponential  decay 
given  by  Eq.  (11a). 

We  may  gain  some  Insight  into  the  general  effect  of  radiation  damping  on 
the  transverse  emittance  by  studying  Eq.  (12)  in  the  case  of  untapered  wlggler 
field.  We  shall  first  consider  the  case  where  Y'  ■  0.  This  could  be  the 
situation  when  the  acceleration  mechanism  is  saturated  by  the  radiation 
damping  and  the  beam  energy  Is  constant.  The  evolution  of  the  beam  radius  Is 
then  given  by  the  appropriate  limit  of  Eq.  (13).  Since  there  is  no  tapering 
of  the  wlggler,  the  solution  is  exact  and  given  by 

a  -  a  1  ♦  sin  2K  z  . 

o  2Kg  B 


The  beam  radius  again  exponentially  decays  with  an  Induced  betatron 
oscillation.  Since  Y  is  constant,  the  damping  rate  v  Is  constant,  and  the 
normalized  emittance  Is  given  by  e^(z)  “  exp(->vz). 

Next,  we  consider  the  situation  when  an  accelerated  beam  Is  cooled  by 
passing  It  through  an  untapered  external  wlggler  field.  Since  the  beam 
decelerates  due  to  the  synchrotron  radiation  damping,  we  have  y'/y  *  -v. 


This  gives  p  -  0  and  since  K_  -  a  k  /(/2Y),  the  betatron  wave  number  Ko  Is  a 

o  w  w  ® 

function  of  z.  The  spatial  dependence  of  Y  can  be  evaluated  using  Y'/Y  -  -v, 

2  2  2  2 

and  Eq.  (13)  reduces  to  a  -  a  (1  ♦  v^z  ),  where  •  Althou^i  the 

O  O  Q  if  W  if  O 

beam  radius  remains  constant  up  to  order  of  0{z  ),  the  normalized  beam 

emlttance  decreases  algebraically,  e  -  e^(0)/(1+v^z). 

n  n  o 

The  relevance  of  the  above  analysis  depends  on  the  magnitude  of  the 
damping  rate  v^.  For  the  following  set  of  accelerator  parameters, ‘ 

Ej^  «  1.5x10^  V/cm.,  -  50  kG.,  -  1  m,  it  is  estimated  that  the  e-fold 

5 

length,  1/Vjj,  could  be  as  short  as  <  600  m  for  Y^  -  10  .  Therefore,  our 
results  show  that  one  can  Improve,  by  Induced  synchrotron  radiation,  the 
quality  of  an  electron  beam  by  passing  it  through  an  external  wiggler  field. 
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IV.  Quantum  Excitation 


An  estimate  for  the  minimum  transverse  normalized  beaun  emlttance  due  to 

quantum  excitation  In  an  IFEL  accelerator  can  be  obtained  from  the  following 

qualitative  treatment.  Similar  arguments  can  be  made  for  electron  beams  In 

storage  rings. The  normalized  transverse  velocity  and  radial 

displacement  of  an  electron  In  a  wlggler  field  are  given  by  fi^-a^/Y,  and 

r^-a^i^/(2»Y).  For  a  fluctuation  6E  In  the  energy  of  the  electron,  the 

corresponding  fluctuations  In  r  and  0  are  6r  -ti6E/E,  and  68  -E6E/E,  where 

n*a^X^/(2iiY)  and  The  Increase  In  normalized  emlttance  due  to  such 

fluctuations  Is]^*^^  Ac^  -  Y[Kg<6r^>  ♦  <6B^>/Kg],  which  for  a  weakly  focusing 

channel,  can  be  approximated  by  Ae^  •  Y<68^>/Kg  -  (YE^/Kg)<6E^>/E^  . 

2 

Due  to  the  discrete  nature  of  the  synchrotron  radiation,  <6E  >  Is  given 
2 

by  N(rtw  )  where  N-Pz/(crt«  )  Is  the  number  of  photons  emitted  In  a  distance  z, 
o  o 

P  Is  the  synchrotron  radiation  power,  and  Is  the  energy  associated  with  a 
quantum  of  synchrotron  radiation.  We  can  therefore  obtain  the  rate  of  change 
of  due  to  quantum  excitation. 


dc  2  Phu 

( _ H)  .  I§ - £ 

^dz  ^Q.E.  Kg  ^g2  • 

However,  with  radiation  damping,  the  total  change  In  is  given  by 


de  dc 

<3r’  * '  '’S  *  <ar>Q.E, 


The  normalized  emlttance,  e^,  reaches  a  minimum,  de^/dz>0,  when  the  two 

2 

effects  are  balanced.  This  gives  minimum  normalized 

emlttance,  where  we  have  used  vc-P/E.  For  synchrotron  radiation, 

■a 

where  p-Y/(a^k^)  Is  the  radius  of  curvature  of  the  electron 
orbit  In  the  wlggler.  The  minimum  transverse  normalized  beam  emlttance  Is 
then  approximately  given  by 
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In  the  case  of  weak  focusing  due  to  wiggler  transverse  gradients. 

Kg  -  a^k^/(v^T),  and  the  nlnlam  nomallzed  emlttance  is 

"n'mln  '  .  (15) 

Using  the  accelerator  parameters  at  the  end  of  section  III,  Eq.  (15) 
gives  the  value  of  the  minimum  normalized  emlttance  to  be  -1.8  cm-rad. 

Such  a  large  value  of  the  minimum  emlttance  Indicates  the  inadequacy  of  the 
weak  focusing  from  the  wiggler  transverse  gradients.  Strong  focusing  from, 
for  example,  a  rotating  quadrupole  field  produced  by  a  pair  of  (or  four) 
helical  current  windings^®* may  be  required.  The  betatron  wavenumber  for 
such  a  focusing  mechanism^'^  is  given  by  Kg  -  je|  OB/3r)/Ym^c  ,  where  3B/3r  is 
the  magnetic  field  gradient  of  the  quadrupole  field  on  axis.  For 

C 

3B/3r-250  G/cm,  a^-600,  A  -10m,  and  T-*lx10  ’  Eq.(l4)  gives  a  minimum 

normalized  emittance  of  e  -0.13  cm-rad.  Another  possible  strong  focusing 

n 

force  could  be  the  electrostatic  radial  electric  field  of  an  ion  coltmm.  Such 

a  column  could  be  created  by  the  ionization  of  the  residual  gas  by  a  low 

energy,  high  current  electron  beam  pulse  preceedlng  the  main  accelerating  beam 

pulse. The  betatron  wavenumber  for  such  a  focusing  mechanism  can  be 

easily  shown  to  be  ,  where  Is  the  Ion  plasma 

frequency  and  (m^^/m^)  is  the  mass  ratio  between  the  ions  and  the  electrons. 

For  n^-10^^/cm^,  a^-600,  A^-10m,  and  T-HxlO^,  Eq.  (1i|)  gives  a  minimum 

normalized  emittance  of  e  -0.04  cm-rad.  An  additional  benefit  of  having  Ion 

n 

focusing  In  the  IFEL  accelerator  Is  that  the  radial  plasma  electron  density 
profile  in  an  ion  column  can  also  be  a  focusing  medium  for  the  laser  beam. 


V,  Nunerlcal  Example 


We  shall  consider  only  resonant  particles  whose  phase  satisfies  the 

2  2 

conditions  d^i/dz  -  0  and  d  i(*/dz  -  0.  The  first  condition  gives 


^  -  l«l/g  B  k 
»  /an  /  ’ 

O 


1/2  4 

y'  -  -  R2  k^"*. 


(16a) 

(16b) 


V 


/k  . 


(I60) 


where  R^  -  /Uje |ELSlnij»jj/(mQC^^I« )  ,  R2  -  je l^k/Om^c^ °)  ,  is  the  resonance 
phase,  the  laser  electric  field  stroigth,  and  k  the  laser  wave  number.  The 
second  condition  together  with  the  pendulun  equation,  Eq.  (6),  provide  the 
spatial  dependences  of  k^  and  B^,, 


3 


k' 


w 


.^sln 


2/^C^  R 


0. 


(17) 


Equation  (17)  shows  that  the  required  tapering  of  the  wiggler  field  may 
be  obtained  by  prescribing  and  a  relationship  between  and  in 
Eq.  (17).  As  an  example,  we  assume  the  tapering  of  the  wiggler  field  to  be 
that  of  a  maxiauB  rate  IFEL  accelerator.*'  For  such  a  case  the  wiggler 
strength  and  the  wiggler  period  are  related  by  the  following  power  law. 


Equation  (17)  nay  then  be  solved  to  give 


B^(0)[1  ♦  R^z]"'^''^, 


(18a) 
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-8/9 


{18b) 


where 


9/2mc'‘ 
lei  /k 


(R^/6R2) B^(0) 


Evaluating  Eqs.  (11)  and  (16a)  with  (16c)  and  (18a, b)  gives  the 
nromalized  transverse  emittance  and  the  resonant  energy  of  the  beam 
as  functions  of  the  propagation  distances. 

For  our  example,  we  ill  consider  the  following  set  of  acceler¬ 


ator  parameters  : 


Ej^  =  1.5  X  10  V/cm,  B^{0)  =  50  kG,  X^{0)  =  100  cm. 


and  X  =  10.6  ^m  with  a  resonance  phase  of  sin  ip-r,  =  0.6.  The  initial 
conditions  are  for  a  matched  beam  with  a  radious  of  1mm  and  a 


normalized  emittance  of  =  0.205  cm-rad,  and  the  required  beam 
injection  energy  is  v52  GeV.  The  beam  is  allowed  to  propagate  for 
1  km  without  depleting  the  laser  radiation.  We  repeated  the  calcula¬ 
tion  by  assuming  there  is  not  radiation  damping  but  with  the  same 
power  law  tapering  of  the  wiggler  field. 

The  results  are  represented  in  Figs.  1,  2,  and  3.  The  open 
squares  denote  the  presence  of  radiation  damping,  while  open 
circles  denote  its  absence.  From  Fig.  1,  we  can  see  that  the 
final  energy  is  not  significantly  reduced  by  the  radiation  damping. 
Figure  2  shows  the  exponential  decay  of  the  normalized  emittance. 

At  the  end  of  the  one-]cilometer  accelerator,  the  normalized 
emittance  is  reduced  to  0.05  cm-rad,  which  is  very  close  to  the 
minimiim  normalized  emittance  of  '^^0.04  cm-rad  at  that  point  if 
ion-column  focusing  is  assumed  in  the  accelerator.  In  Fig.  3,  the 
appropriate  tapering  of  and  for  the  two  cases  is  shown. 
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Conclusion 


We  have  studied  the  evolution  of  transverse  emlttance  and  the  beam  radius 
due  to  the  radiation  damping  effect  In  an  IFEL  accelerator.  We  derived  the 
beam  envelope  equation,  Eq.  (10),  which  Includes  the  effects  of  radiation 
damping,  and  have  demonstrated  that  the  normalized  transverse  emlttance 
decreases  exponentially  with  a  damping  rate  given  by  the  radiation  damping 
coefficient  v  until  It  reaches  a  minimum  value  due  to  quantum  excitation  . 

The  beam  envelope  equation  was  solved  analytically  for  a  slowly-varying 
wlggler  field.  We  have  derived  an  expression  for  the  minimum  normalized 
emlttance  in  the  IFEL  accelerator  and  showed  that  strong  focusing  Is  essential 
In  reducing  this  minimum  emlttance  due  to  quantum  excitation.  We  have  shown 
that  radiation  damping  can  play  an  Important  role  In  Improving  beam  quality 
without  a  significant  sacrifice  in  beam  energy. 
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Figure  Captions 


Fig.  1  Evolution  of  beam  energy  In  an  IFEL  accelerator  with  and  without 
radiation  damping. 

2 

Fig.  2  Exponential  decay  of  normalized  beam  emlttance,  e 

n 

Fig.  3  Spatial  tapering  of  wlggler  period  and  field  with  and  without 
radiation  damping. 
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Appendix  II: 

Radiation  focusing  and  guiding  with  applications 
to  the  Free  Electron  Laser 


10 


In  the  one-diaensional  analysis  of  the  free  electron  laser  (PEL)  the 

radiation  field,  viggler  field  and  electron  beaa  resonantly  couple  so  as  to 

1-3 

Bodify  the  longitudinal  wave  number  of  the  radiation  field.  This 

resonant  interaction  can  lead  to  focusing  of  the  radiation  beam.  This 

phenomena  vas  first  analyzed  for  the  lov  gain  PEL  with  transverse  effects^ 

vhere  it  vas  shown  that  the  diffractive  spreading  of  the  radiation  beam 

could  be  overcome  by  a  focusing  effect  arising  from  the  modified  index  of 

refraction.  This  radiation  focusing  phenomena  has  been  shown  to  play  a 

central  role  in  the  practical  utilization  of  the  PEL,^  since,  in  many 

proposed  experiments  the  radiation  beam  will  not  be  confined  or  guided  by  a 

waveguide  structure.  Recently  optical  guiding  in  PELs  has  been  studied  in 

6-9 

the  small  signal,  exponential  growth  regime,  for  the  asymptotic  behavior 
of  the  radiation  beam. 

In  this  letter,  we  present  a  general,  self-consistent,  fully 
nonlinear,  modal  representation  formalism  which  we  apply  to  the  phenomena 
of  radiation  focusing  and  guiding  in  FELs.  The  novel  aspect  of  our  modal 
expansion  is  that  the  characteristics  of  the  modes  are  governed  by  the 
driving  current  density,  as  opposed  to  a  heuristic  numerical  approach, 
and  hence  it  is  called  the  "source  dependent  expansion"  (SDE).  Instead  of 
using  the  usual  modal  expansion  consisting  of  vacuum  Laguerre-Gaussian 
functions^^  we  incorporate  the  source  function  (driving  current)  self- 
consistently  into  the  functional  dependence  of,  i)  the  radiation  waist,  ii) 
the  radiation  wave  front  curvature,  as  well  as  iii)  the  radiation  complex 
amplitude.  Because  of  the  source  dependent  nature  of  our  modal  expansion, 
the  fundamental  mode  remains  dominant  throughtout  Xhe  evolution  of  the 
radiation  field.  This  approach,  which  can  be  applied  to  a  wide  range  of 
problems,  lends  itself  to  fast  and  accurate  numerical  solutions  as  well  as 
to  a  better  analytical  description  of  the  PEL  focusing  and  guiding  problem. 


1 


Using  the  SDE  approach  in  numerical  siaulations  of  the  PEL,  one  can 
efficiently  incorporate  siaultaneously  the  effects  of  electron  beaa 
eaittance,  energy  spread,  viggler  gradients,  sideband  frequencies,  etc. 

An  envelope  equation  for  the  radiation  is  derived  vhich  describes  the 
transient  as  veil  as  asyaptotic  behavior  of  the  radiation  beaa.  The 
effects  on  the  radiation  beaa  of  a  transversely  displaced  electron  beam  as 
veil  as  a  longitudinally  aodulated  electron  beam  have  also  been  considered. 

In  our  Bodel  the  vector  potential  of  the  radiation  field  is 
^(r,e,z,t)  «  (l/2)A(r,e,z)exp(i(«z/c-wt))e^  -f  c.c.,  vhere  A(r,e,z)  is  the 
complex  amplitude,  u  is  the  frequency  and  c-c.  denotes  the  complex 
conjugate.  The  radiation  field  satisfies  the  reduced  vave  equation, 

F  If  (4)  *  T  “  S  fejaCr,®,.)  -  S(r,e,a),  (1) 

w  r  do  / 

vhere  a(r,0,z)  =  |e  |A/m^c  =  lajexpfi^)  is  the  normalized  complex  radiation 
field  amplitude  and  ve  have  assumed  that  a  ^3a/3z  <<  (o/c.  The  source 
function,  S,  has  the  general  form 

S(r,e,z)  =  (fci/c)^(l-n^(r,e,z,a))a(r,e,z),  (2) 

vhere  n(r,6,z,a)  is  the  complex  index  of  refraction. 

Ve  choose  the  folloving  representation  for  a(r,0,z)  in  terms  of 
associated  Laguerre  polynomials, 

a(r,0,z)  -  X;  E  p(e,z)DP(r),  (3) 

m  p 

vhere  m  and  p  =  0,1, 2,---, 


C  „(0,z)  »  a  (z)cos  p0  +  b  (2)sin  p0, 
ID  y  p  ID  y  p  ID  f  P 


(4a) 

(4b) 


2 


In  Eqs.  (4a, b),  a  (z)  and  b  (z)  are  complex,  r  (z)  is  the  radiation 
a,  p  P  s 

spot  size,  a(z)  is  related  to  the  curvature  of  the  wavefront  and  is  the 
associated  Laj^erre  polynomial.  The  z  dependence  of  these  parameters  will 
be  determined  by  the  source  function  in  Eq.  (1). 

Substituting  (3)  into  (1)  and  using  the  orthogonality  properties  of 
L^,  cos  p6,  and  sin  p6,  ve  obtain. 


where 

A  (z)  «=  r  /r  +  i(2in+p+l)  f{l  +  a^lc/ur^  -  ax  /t  +  o  /2I,  (6a) 

ni,p  ss  \  s  ss  J 

B(z)  =  -  ^ar^/r^  +  (1  -  a^)c/<i)r^  -  a  /2j  -  i  ^r^/r^  -  2ac/(»)r^j,  (6b) 

the  prime  denotes  3/3z,  *  denotes  the  complex  conjugate  and 


2n  ® 

•  2l;(4)T  S(t.e,r)(DP(t))* 

o  o 

(6c, d) 

2  2 

where  C  «  2r  /r  .  The  equations  for  a  and  b  in  (5)  are 
S  D I  P  ^  >  P 

underdetermined,  since  the  function  B(z)  can  be  shown  to  be  arbitrary.  If 
we  choose  B(z)  =  0,  for  example,  we  would  in  effect  be  expanding  the 
radiation  field  in  the  conventional  vacuum  Laguer re-Gauss ian  modes. For 
a  source  free  medium,  B  =  0  would  be  the  most  appropriate  choice.  In  the 
presence  of  a  source  term  a  more  appropriate  choice  for  B(z)  can  be  found. 
This  is  accomplished  by  considering  the  case  where  the  radiation  beam  at 


(1  +  6  _,)”  cos  p6 

p,0‘'  ^ 

sin  p0 


3 


z  -  0  has  a  Gaussian  radial  profile  sjnunetric  about  the  z-axis.  Let  us 

further  assuae  that  for  z  >  0  the  radiation  beaa  profile  remains 

approximately  Gaussian  vith  a  nearly  circular  cross  section.  In  this  case, 

ve  expect  the  magnitude  of  the  coefficients,  a  (z)  and  b  (z)  to  become 

m,p  m,p 

progressively  smaller  as  a  and  p  take  on  larger  values.  A  good 
approximation  to  the  radiation  beam  is  then  given  by  the  lovest  order  mode, 
^0  0^^^'  ve  find  that  only  the  m  -  0,1  and  p  -  0  equations  are 

relevant  and  they  are  (3/3z  +  Aq  q)*q  q  *  -  q  aJ'd  q  ■  ®®q  q*  The 
second  equation  provides  us  vith  a  specific  expression  for  B(z)  in  terms  of 
one  of  the  moments,  q,  of  the  source  term.  The  choice  of  B(z)  «= 

Fj  q(z)  is  source  dependent  and  when  substituted  into  (6b)  yields 

first  order  coupled  differential  equations  for  the  parameters,  r^  and  a,  of 
the  Laguerre-Gaussian  expansion  in  (3)  and  (4a, b).  The  set  of  equations  in 
(5a, b)  may  nov  be  solved  self-consistently  for  the  modal  coefficients  a 

m,  p 

and  b 

m,p 

Ve  first  consider  the  dynamics  of  an  axially  symmetric  radiation  field 
in  the  FEL.  The  appropriate  index  of  ref raction^’ ^ for  a  Gaussian 
beam  density  profile  is 


n(r,z,a) 


1 


1  ^ 

^  2  ^2  V  /  |a(r,2) 


(7) 


2  2  2  2  2 

where  <i)^(r,z)  =  ‘^Q(r^Q/rjj(z) )  exp(-r  /z^{z)),  r^(z)  is  the  electron  beam 
radius,  r^^^  =  is  the  initial  beam  plasma  frequency  on  axis, 

a^  =  |e|B^/k^m^c  is  the  normalized  viggler  amplitude,  y  is  the  electron's 
Lorentz  factor,  i*-  is  the  electron's  phase  in  the  ponderomotive  wave 
potential  and  <  >  denotes  the  ensemble  average  over  all  electrons.  Vith 

the  assumption  that  in  the  source  function  the  complex  radiation  amplitude 


4 


can  be  approximated  by  the  lowest  order  mode,  ve  find  that  (2)  can  be 
written  as 


S(£.,z)  .  -4N)(a^/r^)  , 


^  V  / 


VI 


(8) 


2  3 

where  v  ■  *  Ijj/17xl0  is  Budker's  constant  and  is  the 

electron  beam  current  in  amperes. 

An  envelope  equation  for  the  radiation  beam  can  be  obtained  using  (8) 
and  (6b), 


r^'  +  K^(z,rjj,r^,  Uq  Q|)r^ 


0, 


(9) 


where 


2  2 
=  (2c/«)^ 


-1  +  C^<sin4'>^  +2C<cos<f'>  +  ((i)/2c)r^  C  <sin<k> 

s 


-4 
s  ' 

(10) 


and  C(z)  «(2v/Y)H(z)ay/ |aQ  q(z)|,  measures  the  coupling  between  the 

2  2  2 

radiation  and  electron  beam,  H(z)  -(1-F)/(1+F)  and  F(z)  -  ^b^*^s 

filling  factor.  The  first  term  on  the  right  hand  side  of  (10)  is  the  usual 

diffraction  term,  the  second  and  third  terms  are  always  focusing  while  the 

last  term  is  usually  a  defocusing  contribution.  In  the  high  gain  trapped 

particle  regime,  <sin<</>  and  <cos'k>  are  approximately  constant,  while 

|a  (z)  I  increases  with  z.  Hence,  K  strongly  depends  on  z  and  a  guided 

beam  (r^  =  0)  cannot  be  maintained.  In  the  low  gain  trapped  particle 

regime  |a  (z) |  increases  slightly  and,  therefore,  a  guided  beam  can  be 
0^0 
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approximately  achieved.  In  either  the  Compton  or  Raman  exponential  gain 
regime,  conditions  for  a  stable  guided  beam  can  be  found. 

The  FEL  parameters  used  in  the  following  illustrations  are  similar  to 
those  used  in  Ref.  14  and  are  given  in  Table  I  where  the  resonant  phase 


5 


2 

approximation  is  used  and  z  =  nr  (0)/X  is  the  Rayleigh  length.  Ue  present 

K  S 

first  a  comparison  between;  a)  the  exact  numerical  solution  of  the  wave 

equation  in  (1),  (using  6^x6^  Fourier  modes),  b)  the  solution  using  a 

vacuum  Laguerre-Gaussian  modal  expansion  (B  =  0,  using  10  modes)  and  c)  the 

solution  from  the  Laguerre-Gaussian  SDE  approach  (B  =  q/^q  Qt  using  10 

modes).  For  an  axially  symmetric  configuration,  we  show  in  Fig.  1  the 

evolution  of  the  radiation  beam  amplitude  on-axis  obtained  from  methods 

(a),  (b)  and  (c).  The  SDE  solution  (c)  is  in  excellent  agreement  with 

solution  (a)  while  solution  (b),  beyond  a  Rayleigh  length,  grossly  deviates 

from  (a)  and  (c).  The  excellent  results  obtained  with  the  SDE  approach  are 

also  reflected  in  the  radiation  amplitude  profile.  Figure  2  shows  the 

evolution  of  the  radiation  beam  radius,  r^,  in  the  linear,  exponential  gain 

regime  of  the  FEL  for  the  parameters  in  Table  I.  Five  transverse  modes 

were  used  in  the  numerical  calculation. 

Ue  now  consider  the  case  where  the  electron  beam  centroid  is  displaced 

transversely  in  the  x  direction.  The  index  of  refraction  in  this  case  is 

2  2 

given  by  (7)  with  w^(r,z)  multiplied  by  (1  +  2(r^Xj^/rj^)cos0)  where  is 

the  displacement  of  the  electron  beam's  centroid  and  |xj^|  <<  r^^.  In  the 
FEL  source  term  we  consider  only  the  lowest  order  symmetric  and  anti¬ 
symmetric  mode  with  respect  to  the  x  axis.  The  centroid  of  the  radiation 
beam  is  given  approximately  by 


x^(z)  = 


/2 


(a 


Oil 


va 


0,0'R 


(11) 


II  '  /  / 

where  x,  is  defined  so  that  |a [  is  proportional  to  exp(-((x  -  x,  )"  +  y  )/r  ) 
and  (  )p  denotes  the  real  part.  Figure  3  shows  the  electron  and  radiation 

A 

beam  centroids,  Xj^  =  x^  ( l-sech(k^z) )  and  x^  for  x^  =  =  .075cm 


and  X  =  2ii/k  =  z„/4  =  2.7m.  In  these  numerical  illustrations,  10  radial 
c  c  R 
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nodes  and  2  angular  modes  were  used.  After  an  initial  transient,  the 

radiation  centroid  is  guided  by  and  oscillates  about  the  electron  beam's 

centroid.  We  have  also  studied  the  situation  vhere  the  electron  beam  centroid 

oscillates  according  to  x.  -  x  sin  k  z  with  x  <<  r,  and  X  -  2ri/k  <  2„. 

be  c  cb  c  c-R 

Because  of  the  high  gain  in  the  radiation  field  the  radiation  centroid 
eventually  follows  the  average  position  of  the  electron  beam's  centroid.  When 
the  electron  beam  centroid  oscixlation  is  due  to  the  viggler  field,  there  is 
no  change  in  the  evolution  of  the  radiation  field. 

Under  certain  conditions  the  electron  beam  envelope  can  be  spatially 
modulated  about  the  z-axis  if  the  weak  focusing  force  due  to  the  viggler 
gradient  is  not  balanced  by  the  defocusing  forces  arising  from  emittance  and 
self  field  effects.  It  can  be  shown  that  the  amplitude  and  waist  of  the 
radiation  field  undergo  a  modulation  similar  to  the  electron  beam  envelope 
modulation. 

We  have  analyzed,  using  the  SDE  formalism,  a  number  of  effects  associated 
with  radiation  focusing  and  guiding  in  the  PEL.  This  approach  can  be  readily 
generalized  to  include  both  spatial  and  temporal  variations  in  the  radiation 
field  in  order  to  study  sideband  generation  and  focusing  effects 
simultaneously  in  the  PEL. 
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Table  I 


Electron  Bean 


Current 

1 

M 

2kA,  (V 

-  0.118) 

Energy 

50  MeV, 

(y  -  100) 

Radius 

-  0.3  cn 

Radiation  Beam 

Wavelength 

X  -  10.6pm 

Input  Power 

P(2.0)  -  230MV,  (|a(0,0)|  -  1.84x10“'^) 

Spot  Size 

r  (0)  «  0.6  cm,  (Zp  «  10.7  m) 

S  K 

Viggler  Field 

Wavelength 

X  =  8  cm 

V 

Wiggler  Strength 

B^  =  2.3  kG,  (a^  -  1.716) 

Resonant  Phase 

»=  0.358  rad 
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Figure  Captions 


Fig.  1  Radiation  amplitude  on  axis,  |a(0,z)|  for  a)  exact  numerical 

solution  (64x64  Fourier  modes),  b)  vacuum  modal  expansion  solution 
(10  modes),  and  c)  SDE  solution  (10  modes)  at  distance  of 
z  «  42-  «=  42.8  m. 

K 

Fig.  2  Evolution  of  the  radiation  beam  radius,  1/e  vidtn.  for  initial 

spot  sizes:  a)  0.35  cm,  b)  0.24  cm,  and  c)  0.15  cm. 

Fig.  3  Electron  and  radiation  beam  centroids,  and  x^  for  a  displaced 
electron  beam,  x.  =  x  (l-sech(k  z))  with  x  r./4  and 

DC  C  CD 

X  =  2ii/k  = 
c  c  R 
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Table  Caption 


Table  1. 


PEL  simulation  parameters 


7X10 


Z/Z 


z  (m) 


/z 


Appendix  III: 

Analysis  of  radiation  focusing  and  steering  in 
the  Free  Electron  Laser  by  use  of  a  source  dependent 
expansion  technique 
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II ,  Foripulatlon  of  the  Source  Dependent  Expansion  (SDE) 


The  radiation  focusing  and  guiding  configuration  for  the  FEL  is  shovn 
in  Fig.  1.  In  our  model  the  vector  potential  of  the  linearly  polarized 
radiation  field  is 


^(r,e,z,  t) 


_  A(r,e,z) 
2 


^i  ( (oz/c-wt ) ' 


c.  c. 


(1) 


where  A(r,e,z)  is  the  complex  radiation  field  amplitude,  «  is  the  frequency 
and  c.c.  denotes  the  complex  conjugate. 

The  wave  equation  is 


1  3_ 

fr  i-1 

2  2 

1  3  3 

2  1 
1  3^ 

L  - 

in 

r  3r 

1  3rJ 

*  1  1  *  2 
r  36  3z 

c  3t  J 

-R  “ 

c 

J  e  , 
x  X 


(2) 


where  J^(r,0,z,t)  is  the  driving  current  density  associated  with  the 
medium.  Substituting  (1)  into  (2)  leads  to  the  following  reduced  wave 
equation, 


i_  i_ 
+  22 
r  36 


2i 


c 


\ 

8_ 

3z 


a(r, e,z) 


j 


S(r , 0, z) , 


(3) 


2 

where  a(r,6,z)  =  |e|A/m^c'  =  jalexpCi^)  is  the  normalized  complex  radiation 
field  amplitude  and  we  have  assumed  that  a(r,6,z)  is  a  slowly  varying 
function  of  z,  i.e.,  a~^3a/3z  <<  w/c.  The  amplitude,  |a(r,6,z)|,  and  phase 
♦(r,0,z)  are  real  functions  expressed  in  terms  of  the  polar  coordinates, 
r,6  and  z.  The  source  function,  S,  has  the  general  form 

2 

S(r,6,z)  =  ^  (l-n'^(r,0,z,a))a(r,6,z),  (4) 

c 


where  n(r,6,z,a)  is  the  index  of  refraction  associated  with  the  medium  and 
is  in  general  complex  and  a  function  of  r,6,z  as  well  as  the  radiation 
field,  a. 


4 


Ve  choose  the  following  representation  for  a(r,e,2)  In  terms  of 
associated  Laguerre  polynomials, 


a(r,e,z)  =  E  E  C  (e,z)DP(r), 

m ,  p  m 


where  m  «  0,1,2,*--,  p  =  0,1,2,'-*, 


C  (6,z)  =  a  (z)cos  p0  +  b  (z)sin  p6, 
m.p  ni,p''  m,p'  ^ 


In  Eqs.  (6a, b),  the  complex  coefficients  a  (z)  and  b  (z)  are  functions 

of  2,  is  the  radiation  spot  size,  a(z)  is  related  to  the  inverse  of 

the  radius  of  curvature  of  the  radiation  beam  (curvature  of  wavefront)  and 

lP  is  the  associated  Laguerre  polynomial.  Solving  for  the  unknown 

quantities  a  ,  b  ,  r  and  a  in  terms  of  the  source  term  S  allows  us  to 
m,p’  m,p’  s 

completely  describe  the  radiation  dynamics.  It  will  be  shown  later  that 
the  representation  in  (5)  is  underspecified,  there  are  more  unknown 
quantities  in  (5)  than  available  equations.  The  additional  degrees  of 
freedom  in  our  representation  allow  us  to  specify  particular  functional 
relationships  for  the  unknown  quantities  r^  and  o  in  such  a  way  that  the 
number  of  terms  (modes  needed  to  accurately  describe  the  radiation  beam)  is 


small . 


To  proceed  with  the  derivation  we  substitute  (5)  into  (3)  and  obtain 


3C  (e,z)/3z  C  (e,z)<^ 
m,p  ’  '  m,p  |3z 


.  ^  fi-  fr 

u)r_  )) 


2w  S(£.,e,2), 
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2  2 

vhere  *  2r  /r  (z).  It  can  be  shown  that  the  second  term  on  the  left  side 
s 

of  (7)  can  be  put  into  the  form 


'  -  i(m+l)B(z)DP^j(0  -  i  (m4p)B*(z)DP_  j(  O  ,  (8) 

vhere 

^  i(2m+p+l)^(l  +  a^)c/wr^  -  <xr^/r^  +  o  /2j,  (9a) 

B(z)  =  -  ^ar^/r^  +  (1  -  a^)c/ur^  -  a  /2j  -  i  ^r^/r^  -  2ac/a)r^j,  (9b) 

*  denotes  the  complex  conjugate  and  the  prime  denotes  a  derivative  with 
respect  to  z,  i.e.,  '  =  3/3z. 

In  obtaining  (8)  the  following  identities  were  used: 

tuP  .  (2.»p*l)uP  -  (.*1)UP^J  -  <n,*p)uP_j, 

2ldU^/di  =  (2ni+p-C)U^  -  2(m+p)uP  and 
m  tn  m- 1 

£.3^uP/3£;^  +  -  2(2m+p+l)juP, 

vhere  U^ff;)  =  ?;)exp(-£,/2) .  Substituting  (8)  into  (7)  and  performing 

the  operation 
2n 

J  (cos  p'G,  sin  p'0)d0/2n, 
o 

on  the  resulting  equation  yields 
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Ik 


3/3z 


m,pj 


m=0 


'a  1 

A  ^  V  1 

ra  'll 

m,p 

b 

-  ‘  “Li  ‘  Cl 

cr 

B 

V  m, pj 

1 

I.  ni,p;j 

zis.  r 

2ria)  J 


2n 

de  S({.,e,z) 


pel 

sin  pe 


(10a, b) 


where  6^  ^  is  the  Kronecker  delta.  Multiplying  (10)  by  (0^)*  and 
integrating  over  (.  from  0  to  «  yields 


f-  +  A  (z) 
3z  m,p 


m,p 

3 

m.P 


-  imB(z) 


m-l,p 


(2) 


lb  ,  (z) 

m-1 ,  p 


-  i(m+p+l)B  (z) 


a  1 

m+1 , p  ' 

|b  ,  (z) 

(  m+l,p' 


=  -  1 


(T  (2)) 

m.P 

G  (z)  ’ 

.  m.P  ) 


(11a, b) 


where 


2n 


^(U6p^Q)"kos  pO) 

sin  p9 


(12a, b) 


In  obtaining  (11)  we  used  the  orthogonality  relation. 


jDP<i;)(DP(t)]*dt.  12^  «, 

O 


The  function  B(z)  is  arbitrary  and  is  not  specified.  The  equations 


for  a  and  b  in  (11)  are  underdetermined,  since  the  function  B(z)  can 
m ,  p  m ,  p 

be  shown  to  be  arbitrary.  If  we  choose  B(z)  =  0,  for  example,  we  would  in 


effect  be  expanding  the  radiation  field  in  the  conventional  vacuum 
Laguerre-Gaussian  modes. Ue  will  show  later  that,  in  general,  expansion 
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in  terms  of  the  vacuum  modes,  B  -  0,  would  require  far  too  many  modes  to 

accurately  describe  the  radiation  beam  over  distances  of  many  Rayleigh 

lengths.  A  more  appropriate  choice  for  8(2)  will  depend  on  the  particular 

problem  under  consideration.  Let  us  consider  one  of  the  most  common 

situations  where  the  radiation  beam  at  2  ■=  0  is  known  and  has  a  Gaussian 

radial  profile  symmetric  about  the  z-axis.  In  this  case  the  complex 

radiation  amplitude  at  z  •=  0,  is  given  by  a(r,0,O)  '  ^Iq  q 
7  7 

(-(1  -  ia(0))r  /r^(0))  and  is  independent  of  0.  Let  us  further  assume  that 
for  2  >  0  the  radiation  beam  profile  remains  approximately  Gaussian  with  a 


nearly  circular  cross  section.  That  is,  the  dominant  part  of  the  source 

S(r,0,2)  has  an  r  and  2  dependence  and  the  0  dependent  part  is  weak.  In 

this  case  ve  expect  the  magnitude  of  the  coefficients,  a  ^(z)  and  b  „(z) 

,  p  m ,  p 

to  become  progressively  smaller  as  m  and  p  take  on  larger  values,  i.e.. 


la  1  »  la  .  I,  la  „  ,1  and  |b„  »  jb  -  I,  |b  J.  The  lowest 

'  m,p'  '  m+l,p'  '  m,p+l'  '  in*P  m+l,p‘  m,p+l' 

order  approximation  to  the  radiation  beam  is  given  by  the  a^  q(z)  mode. 


Hence,  if  the  a^  g  mode  gives  a  rough  approximation  to  the  radiation  field 
ve  may  solve  for  a^  q(2)>  «(z)  using  (11a).  From  (11a)  ve  find 

that  only  the  m  =  0,1  and  p  =  0  equations  are  relevant  and  yield 


(3/az  +  Ao,0^^0,0  "  "  ^^0,0’ 


(13a) 


®®0,0  "  ^1,0' 


(13b) 


Ve  nov  have  a  specific  expression  for  B(z),  from  (13b),  in  terms  of  one 
of  the  moments,  F,  of  the  source  term.  Substituting  (9b)  into 

1  f  u 

B(z)  =  F^  g(z)/ag  q(z)  yields  the  folloving  first  order  coupled 
differential  equations  for  r^  and  a 

r'  -  2ccx/cor^  -  -  ^5(^1  ,o/^0,0>I  ’ 
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-  2(1  •  -  2((''i,0'“0,0)r  -  “<fl,o(“0,0)l]’  (’‘"I 

vhere  (  )d  t  denotes  the  real  and  imaginary  part  of  the  enclosed  function. 

Since  r  (z)  and  0(2)  are  now  known  from  (14a, b)  we  may  solve  for  A  (z) 

S  ®  I  P 

using  (9a), 


A  (2)  =  2ca/«r  -(P,  ^/a^. 
m,p  s  1,0  0,0  I 


+  1 


i(2m+p>l)(2c/«r^  4  (Fi,o^^O,0^r)  ’ 


(15) 


Using  B(2)  =  Fj  Q(z)/aQ  q(z)  and  the  resulting  equations  for  and  a  in 

(14)  allows  us  to  solve  for  a  and  b  in  (11a, b). 

m,p  m,p 

It  is  useful  at  this  point  to  consider  the  simple  case  of  propagation 

of  a  radiation  beam  in  vacuum  (no  source  term).  To  illustrate  this  well- 

known  limit  we  evaluate  a  ,  b  ,  r  and  a  in  the  source-free  case, 

m,p’  m,p’  s 

' '  2  3 

F  =  G  =  B  =  0.  Equations  (14a, b)  become  r^  =  (2c/w)  /r^  and 
ITllPniyP  3  s 

a  =  (w/2c)r^r^  and  have  the  solutions 


2  2  1/2 

r^(z)  =  r^(0)(l  >  zV2^)"'^ 


(16a) 


a(z)  =  z/Zj^, 


(16b) 


vhere  r  (0)  is  the  minimum  radiation  spot  size  at  z  =  0,  z^  =  (u/2c)r  (0)  = 

S  K  S 

2 

nr^(0)/X  is  the  Rayleigh  length  and  X  =  2nc/«  is  the  wavelength.  From 

2 

(10a)  we  find  that  A  (z)  =  2(o(z)  +  i(2m+p+l))cwr  (z)  which  allows  us  to 

m,p  s 

solve  for  a  and  b  using  (11) 
m,  p  m ,  p 


a  ( z )  ^ 
m,p' 

,b  (z) 

V  m.p' 


a  (0)1  . 

b'"’^0) 

1  m ,  p  2 


-i (2m+p+l ) tan  ^(z/z^) 


(17) 


Equations  (16a, b)  and  (17)  together  with  the  representation  in  (5),  (6a, b) 
is  in  agreement  with  the  conventional  vacuum  Gauss ian-Laguerre  form. 
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Ill .  Radiation  Focusing  and  Steering  in  FELs 


A.  Radiation  Beam  Envelope  Equation 

Ve  first  consider  the  dynamics  of  an  axially  symmetric  radiation 

field  in  the  PEL.  For  a  linearly  polarized  viggler  field  and  axially 

symmetric  electron  beam  having  a  Gaussian  density  profile,  the  appropriate 

4  5  9  12 

index  of  refraction  for  the  PEL  mechanism  '  '  ’  is 


n(r,z,a) 


1 


1  \ 

*  2  ^2  \t  /  |a(r,2)|  ’ 


(18) 


2  2  2  2  2 

where  w^(r,z)  =  “jjq( )  exp{-r  /r^^(2)),  r^^Cz)  is  the  electron  beam 

2  1/2 

radius,  r^^  =  rj^(O),  =  (4n|e|  n^Q/m^)  is  the  initial  beam  plasma 

frequency  on  axis,  n^^^  is  the  initial  beam  density  on  axis, 
a^  =  |e |B^/k^mQC  is  the  normalized  viggler  amplitude,  B^  is  the  viggler 
magnetic  field  strength,  k  is  the  viggler  wave  number,  y  is  the  electron's 
Lorentz  factor,  »</  is  the  electron's  phase  in  the  ponderomoti ve  wave 
potential  and  <  >  denotes  the  ensemble  average  over  all  electrons. 

Substituting  (18)  into  (4)  and  noting  that  ll-n|<<l,  gives  the  PEL  source 
function 


S(r,z) 


-  a^(r,z) 
2 


a(r,z) 

a(r,z) 


(19) 


Since  the  electron  beam  radius,  r^,  may  not  be  matched  with  respect  to  the 

focusing  fields  (viggler  gradients)  and  defocusing  effects  (beam  emittance) 

ve  allow  r^  to  be  a  function  of  z  (this  case  is  considered  in  Sec.  IIIC). 

To  proceed  with  the  analysis  ve  assume  that  in  the  source  function  the 

complex  radiation  field  amplitude  in  (5)  can  be  approximated  by  the  lowest 

2  2 

order  mode,  a^  Q(z)exp(-(l-ia)r  With  this  assumption  the  source 

function  can  be  written  as 


10 


S(C.z)  -  -4v(a^/rJ) 


^0,0  /e' 

''  V 


(20) 


2  3 

where  v  -  '  Ij^/17xl0  is  Budker's  constant  and  is  the 

electron  beam  current  in  amperes.  The  moments  of  the  source  function, 
F  .(z),  are  given  by  (12a) 

ID  f  p 


F  _  =  -  4  -  v(a  /r?)  , 
m,0  <*)  V  b  |a 


.  2  ,  2  , .  m 

0^  /.-A  w> 

^  V  /  (r^/r^4 


1) 


m+1 


,  (21) 


where  we  have  assumed  ij*  to  be  constant  across  the  electron  beam.  Since  we 

are  considering  an  axially  symmetric  electron  beam  and  radiation  field  we 

note  that  a  =  F  =  G  =0  for  p  >  0.  Substituting  (21)  into  (14a, b) 
m,pm,pm,p 

and  (15)  yields 


r  r  -  2c(x/(ti  =  -2-  C(z)  <sin  «k>, 
s  s  0)  '  ' 


(22a) 


u 


r^a  -  2(1  +  o^)c/«  =  -4^  C(z)  ^<cos  +  a  <sin  (22b) 


Am  q(z)=  ^ — l^a  +  i(2m+l)  -  C(z)  ^<sin  +  i  ( 2m+ 1 )  <cos\|/>  j 


r  o) 
s 


(22c) 


where  C(z)  =(2v/Y)H(z)a^/ ja^  g(2)|*  H(z)  =(1-F)/(1  +  F)^  and  F(z)  = 

the  filling  factor.  The  function  C(z)  measures  the  coupling  between  the 

radiation  and  electron  beam  and  decreases  as  the  radiation  grows. 

Equations  (22a)  and  (22b)  can  be  combined  to  give  the  following 
envelope  equation  for  the  radiation  beam 


K  (2,r^,r^,  |ag_o|)r^  =  0, 


(23a) 


where  the  initial  condition  on  r'  is  found  from  (22a)  and 

s 


11 


2  2 
K  .  (2c/u) 


-1  4  C^<sin«^/>^  ■f2C<cos(t'>  4  («/2c)r^  C'<sin«^> 


s 

(23b) 


The  first  term  on  the  right-hand  side  of  (23b)  is  defocusing  and 
corresponds  to  the  usual  diffraction  expansion,  the  second  and  third  terms 
are  always  focusing  while  the  last  term  Is  a  defocusing  contribution. 

B.  Radiation  Focusing 

2 

Focusing  occurs  when  K  >  0.  In  the  high  gain  trapped  particle 

regime,  the  condition  for  a  perfectly  guided  beam  (K  =  0)  cannot  be 

2 

maintained  since  K  decreases  as  the  radiation  grows.  In  the  small  signal, 

exponential  gain  regime  the  quantities  <sini|(>  and  <cos'i'>  may  be  calculated 

from  the  linearized  orbit  equations.  The  envelope  equation  may  then  be 

solved  to  determine  r  as  a  function  of  distance  along  the  viggler.  One 

s 

finds  that  in  this  regime,  conditions  for  a  perfectly  guided  radiation  beam 
can  be  achieved. 

Using  (11a)  or  (13a)  ve  find  that  the  magnitude  of  aQ  q(z)  evolves 
according  to 


0/3z  4  (A 


0,0 


^  ^0,0^^ 


0,0 


"  "^^^0,0  ^0,0 


^0,0®0,0^‘ 


(24) 


Substituting  (21)  and  (22c)  with  m  =  0  into  (24)  and  using  (22a)  yields 


3-  ,  I  h  .  ^  i 
9z  ^^s  *^0,0  I  u  Y  ^w 


,  2  2, 


.2  . 


where  (J^^laQ  q|)  is  proportional  to  the 
(  |aQ  q( z )  |r^/X)^ (Uat ts ] .  Equation  (25) 
(23a)  and  show  that  the  maximum  rate  of 


r  =  r,  . 
s  b 


radiation  power, 
should  be  solved 
increase  in  power 


(25) 

P(z)  =  2.15x10^^ 
together  with 
occurs  when 


12 


C.  Radiation  Steering  in  the  FEL 

In  the  FEL  the  centroid  of  the  electron  beam  may  be  displaced  off- 

axis  by  a  misalignment,  a  redirection  of  the  beam  or  because  of  the 

oscillations  in  the  viggler  field.  To  determine  the  degree  to  vhich  the 

radiation  beam  will  follov  or  be  steered  by  the  electron  beam,  ve  consider 

the  case  where  the  electron  beam  centroid  is  displaced  transversely  in  the 

X  direction.  The  index  of  refraction  in  this  case  is  given  by  (18)  with 
2 

(*J^(r,2)  replaced  by 

-r^/r^ 

w^(r,e,z)  =  ®  ^  (2rXj^(z)/r^)cosej ,  (26) 


where  Xj^(z)  is  the  displacement  of  the  electron  beam's  centroid  and 

|x^|<<  r^^.  In  the  source  term,  given  by  (19),  ve  consider  only  the  lowest 

order  symmetric  and  anti-symmetric  mode  with  respect  to  the  x  axis, 

1/2 

a(r,e,z)  -  (Bq  0  ^0  1  ^  cose)exp(-(  1-ia)  1/2) .  With  this  assumption 

the  moments  of  the  source  function,  F  (z),  for  p  =  0,1  are, 

m,  p 


F  =  -8  -  v(a 
m ,  p  u  V  b  \y 


0,0 


"0,0' 


e(z)  +  i 

(H] 

< 

(rhrl  -1)' 


„  „  m+p+1  ’ 

(27) 


(rVr^  ^1) 


1/2  2 

where  e(z)  =  2  x. (z)r  (z)/r,  and  G  =0.  For  small  displacements  of 

b  s  b  m,p 

the  electron  beam  centroid  it  is  easy  to  show  that  the  centroid  of  the 
radiation  beam  is  given  by 


x^Cz)  = 


i'.-(z)  fa 


/2 


0,0'R 


where  x  is  defined  so  that  |a |  is  proportional  to  exp(-((x 


(28) 

2 


2x/  2, 

y  )/r^). 
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D.  Effect  of  a  Modulated  Electron  Beam 

The  electron  beam  envelope  in  the  PEL  can  undergo  modulations. 

The  modulation  is  symmetric  about  the  z-axis  and  can  be  caused  by  improper 
values  for  the  beam  emittance,  radius  and/or  current  injected  into  the 
viggler  region.  For  small  perturbations  about  the  matched  beam  radius, 
tbo*  find  from  the  electron  beam  envelope  equation  that  r^^fz)  »  r^^^ 

(1  +  Asin(K„z))  vhere  r,„  -  (2e  /a  k.  ,  Ko  «  a  k  /-iZy  is  the  betatron 

D  uU  n  V  V  D  V  V 

wave  number,  due  to  the  veak  focusing  effect  of  viggler  gradients,  is 

2 

the  normalized  emittance,  a  =  le |B  /(k  m  c  )  and  A  <<  1.  The  modulation 
of  the  electron  beam  envelope  may  be  included  in  the  source  term,  Eq.  (19), 
through  the  electron  beam  plasma  frequency,  a5j^(r,z).  The  effect  of  a 
modulated  electron  beam  on  the  radiation  beam  is  illustrated  in  the  next 
section. 

In  cases  vhere  the  electron  beam  centroid  or  envelope  is  displaced 

or  modulated  vith  a  spatial  period  close  to  the  viggler  period,  it  becomes 

necessary  to  include  in  the  source  function,  Eq,  (19),  the  rapidly  varying 

part  of  the  phase  This  rapidly  oscillating  contribution  to  the  phase, 

2  2 

(a  /(4  *  2a  ))sin  2k  z,  arises  from  the  linearly  polarized  viggler  field. 

V  V  V 
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IV.  Numerical  Results 

In  this  section  ve  apply  the  SDE  formulation,  given  by  (11)  together 

with  (12)  to  the  FEL.  Using  the  source  term  given  in  (4)  and  (18)  ve  first 

present  a  comparison  between;  a)  the  exact  numerical  solution  of  the  wave 

equation  in  (4),  (using  64x64  Fourier  modes),  b)  the  solution  obtained 

using  a  vacuum  Laguerre  modal  expansion  (B=0,  using  10  modes)  and  c)  the 

solution  obtained  from  the  Laguerre  SDE  approach  (B  =  F^  q/^Lq  q»  using  10 

modes).  The  FEL  parameters  used  in  these  illustrations  are  similar  to 

those  used  in  Ref.  14  and  are  given  in  Table  I  where  the  resonant  phase 

approximation,  <exp(-i4>)>  =  expf-i^i^),  is  used  for  demonstration  purposes. 

Propagation  distances  are  measured  in  terms  of  the  Rayleigh  length,  = 

K 

2 

nr  (0)/X  where  X  is  the  wavelength  and  r  (0)  is  the  minimum  spot  size. 

w  S 

For  an  axially  symmetric  configuration,  Fig.  2  shows  the  radiation 
magnitude,  |a(r,z)  |  as  a  function  of  r  at  four  Rayleigh  lengths 
(z  =  4Zj^  =  42.8m)  for  the  (a),  (b)  and  (c)  methods  of  solution.  The  SDE 
solution  (c)  shows  excellent  agreement  with  solution  (a),  while  solution 

(b)  is  in  poor  agreement.  To  continue  this  comparison  ve  show  in  Fig.  3 
the  evolution  of  the  radiation  beam  amplitude  on-axis  obtained  from  methods 
(a),  (b)  and  (c)  as  a  function  of  propagation  distance.  The  SDE  solution 

(c)  is  again  in  good  agreement  with  solution  (a)  where  as  solution  (b), 
beyond  a  Rayleigh  length,  grossly  deviates  from  (a)  and  (c).  The  results 
in  Figs.  2  and  3  clearly  show  the  improved  accuracy  of  the  SDE  approach 
over  the  conventional  vacuum  expansion  method.  As  an  example  of  radiation 
focusing,  for  an  FEL  in  the  small  signal,  exponential  gain  regime,  the 
radiation  beam  radius  is  found  to  asymptotically  approach  a  matched 
perfectly  guided  value  as  shown  in  Fig.  4.  For  this  example,  the 
parameters  of  Table  I  were  used  and  five  modes  employed. 
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Ve  nov  use  the  SDE  method  to  illustrate  the  steering  of  the  radiation 

beam  when  the  electron  beam  is  displaced  off-axis.  In  these  numerical 

illustrations,  10  radial  modes  (m  =  0,...,9)  and  2  angular  modes 

(p  =  0,1)  were  used.  In  the  first  example,  the  electron  beam  centroid  is 

displaced  off-axis  according  to  x^  =  x^  ( l-sech(k^z ) ) .  Figure  5  shovs  the 

electron  and  radiation  beam  centroids  x,  and  x,  for  x  «  -  0.075cm  and 

b  L  c  b 

X  =  2n/K  =  42„  =  42.8m.  The  radiation  centroid  follows  and  oscillates 
c  c  R 

about  the  electron  beam's  centroid.  Figure  6  shows  the  radiation  profile 
at  twelve  Rayleigh  lengths  (z  =  12  Zj^).  The  asymmetry  of  the  radiation 
profile  is  apparent.  Figure  7  shows  another  illustration  of  steering  where 
the  electron  beam  is  displaced  more  abruptly,  with  X  =  z_/4  =  2.7m.  After 
an  initial  transient,  the  radiation  centroid  is  again  steered  by  and 
oscillates  about  the  electron  beam's  centroid.  In  the  next  example  we  take 

the  electron  beam  centroid  to  be  oscillating  about  the  z-axis,  x^  =  x^sin 

k^z,  with  amplitude  x^  =  r^j/4  and  period  =  10.7m.  Figure  8  shows 

the  electron  and  radiation  centroids  Xj^  and  Xj^,  as  a  function  of  z/z^. 

Because  of  the  high  gain  in  the  radiation  field,  the  radiation  centroid 
eventually  follows  the  average  position  of  the  electron  beam's  centroid. 
Figure  9  shovs  the  distortion  of  the  radiation  profile  due  to  the 
oscillating  electron  beam  at  twelve  Rayleigh  lengths  (z  =  12Zj^).  In  the 
case  where  the  electron  beam  centroid  oscillation  is  due  to  the  viggler 
field,  X  =  a  /rk  and  k  =  2il/X  ,  no  noticeable  change  in  the  evolution  of 
the  radiation  field  (compared  to  the  case  for  x^  =  0)  is  observed. 

The  last  illustration  is  for  the  case  where  the  electron  beam  envelope 
is  spatially  modulated.  Using  the  parameters  in  Table  I  ve  find  that 

=  0.06  cm-rad  and  Xg  =  2n/Kg  =  4.66m.  Figure  10  shovs  the  amplitude  of 
the  radiation  field  on-axis  as  a  function  of  propagation  distance  when  the 


electron  beam  envelope  is  not  matched,  r,  =  r,„  (1  +  Asin(K_z)),  where 


*^b0  ~  0-3cm  and  A  =  0.1. 
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V.  Conclusion 


In  this  paper  a  technique  for  solving  the  three-dimensional  vave 
equation  vith  a  driving  current  density  has  been  developed.  Using  this 
source  dependent  expansion  technique,  a  number  of  effects  associated  vith 
radiation  focusing  and  steering  in  the  FEL  have  been  illustrated.  The 
formalism  Is  used  to  derive  a  general  envelope  equation  for  the  radiation 
beam.  Using  the  envelope  equation,  ve  find  that  it  is  possible  to  have  a 
stable  guided  optical  beam  in  the  exponer.tial  gain  (small  signal)  regime 
but  not  in  the  high  gain  trapped  particle  regime.  Ve  also  considered  the 
effects  on  the  radiation  beam  vhen  the  electron  beam  centroid  is 
transversely  displaced  and  vhen  the  electron  beam  envelope  is  modulated. 

The  source  dependent  expansion  approach  lends  itself  to  fast  and  accurate 
numerical  solutions  as  well  as  to  a  better  analytical  description  of 
focusing  and  steering  in  the  FEL.  Ve  conclude  by  noting  that  this  approach 
can  be  readily  generalized  to  include  both  spatial  and  temporal  variations 
in  the  radiation  field  in  order  to  study  sideband  generation  and  focusing 
effects  simultaneously  in  the  FEL. 
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Table  I 


Electron  Beam 


Current 

=  2kA,  (V 

-  0.118) 

Energy 

=  50  MeV, 

(y  =  100) 

Radius 

Tbo  =  0.3  cm 

Radiation  Beam 

Wavelength 

X  =  10. bum 

Input  Power 

P(2=0)  =  230MV,  (|a(0,0)|  =  1.84x10"^) 

Spot  Size 

r  (0)  =  0.6  cm,  (z^  =  10.7  m) 

S  K 

Viggler  Field 

Wavelength 

X  =  8  cm 

V 

Wiggler  Strength 

B  =  2.3  kG,  (a  =  1.716) 

V  V 

Resonant  Phase 

=  0.358  rad 
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Figure  Captions 


Pig.  1  Schematic  of  radiation  focusing  and  guiding  in  an  FEL. 

Fig.  2  Radiation  amplitude  profile,  |a(r,z)l  for;  a)  exact  numerical 

solution  (64x64  Fourier  modes),  b)  vacuum  modal  expansion  solution 
(10  modes),  and  c)  SDE  solution  (10  modes)  at  a  distance  of  z  -  4zj^ 
s  42.8  m ■ 


Fig- 


Fig- 


Fig. 


Fig. 

Fig. 


Fig. 

Fig. 

Fig. 


3  Radiation  amplitude  on  axis,  (a(0,z)  |  for;  a)  exact  numerical 
solution  (64x64  Fourier  modes),  b)  vacuum  modal  expansion  solution 
(10  modes),  and  c)  SDE  solution  (10  modes). 

4  Spatial  evolution  of  the  radiation  spot  size  in  the  exponential 
gain  regime  for  initial  spot  sizes;  a)  0.35  cm,  b)  0.24cm,  and 
c)  0.15  cm. 

5  Electron  and  radiation  beam  centroids,  x^^  and  for  a  displaced 
electron  beam,  x.  =  x  (l-sech(k  z))  with  x^  =  r, /4  and 

DC  C  CD 

X  =  2n/k  =  4Zo. 
c  c  R 

6  Radiation  amplitude  profile  at  z  =  12zj^  for  a  displaced  electron 
.  beam,  x^  =  x^( l-sech(k^z) )  vith  x^  =  rj^/4  and  =  2ii/k^  =  4Zj^. 

7  Electron  and  radiation  beam  centroids,  Xj^  and  Xj^  for  a  displaced 
electron  beam,  x,  =  x  (l-sech(k  z))  vith  x  =  r,/4  and 

DC  C  CD 

X  =  2n/k  =  Zi,/4. 
c  c  R 

8  Electron  and  radiation  beam  centroids,  x^  and  x^^  for  an  oscillating 

electron  beam,  x.  =  x  sin  k  z  vith  x  =  r,/4  and  X^  =  2n/k^  =  Zp. 

bcc  cb  c  cK 

9  Radiation  amplitude  profile  at  z  =  12  z^  for  an  oscillating 

electron  beam,  x.  =  x  sin  k  z  vith  x  =  r,/4  and  X  =  2n/k  =  z_. 

bcc  cb  c  cK 

10  Radiation  amplitude  on  axis,  )a(0,z)|  for  a  modulated  electron 
beam,  r^^  =  Tj^qCI  *  ^  sin(KgZ))  vith  r^^  =  O.lcm,  A  =  0.1, 

Xg  =  2R/Kg  =  4.66  m. 
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FEL  Interaction  Region 
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II .  Formulation  of  the  Source  Dependent  Expansion  (SDE) 

The  radiation  focusing  and  guiding  configuration  for  the  FEL  is  shown 
in  Fig.  2.  In  our  model  the  vector  potential  of  the  linearly  polarized 
radiation  field  is  A_(r,e,2,t)  =  ( 1/2 )A( r , 6, z )exp(  i  ( wz/c-wt ) )e  +  c.c., 

•“K  X 

where  A(r,0,z)  is  the  complex  radiation  amplitude,  w  is  the  frequency 

and  c.c.  denotes  the  complex  conjugate. 

2-2  2  2 

The  radiation  field  satisfies  the  wave  equation  (V  -c  8  /8t  )Aj^  = 
-4nc~^J^e^,  where  J^(r,G,z,t)  is  the  driving  current  density  associated  with 
the  medium.  Substituting  A„  into  the  wave  equation  leads  to  the  following 
reduced  wave  equation, 


1  ^ 
r  3r 


u 

c 


dz 


a(r, e,z) 


S(r,e,2), 


(1) 


2 

where  a(r,e,z)  =  |e|A/m^c  =  jalexpCiij))  is  the  normalized  complex  radiation 
field  amplitude  and  we  have  assumed  that  a(r,0,z)  is  a  slowly  varying 
function  of  z,  i.e.,  a”^3a/3z  <<  w/c.  The  source  function,  S,  has  the 
general  form 


S(r,0,z) 


~  (l-n‘^(r,0,z,a))a(r,0,z), 
c 


(2) 


where  n(r,0,z,a)  is  the  index  of  refraction  associated  with  the  medium  and 
is  in  general  complex. 


Ue  choose  the  following  representation  for  a(r,0,z)  in  terms  of 
associated  Laguerre  polynomials, 

a(r,0,z)  =  Z  E 
m  p 


3 


where  m  =  0,1,2, •••,  p  =  0,1,2,‘*', 


C  (6,z)  =  a  (z)cos  p0  4  b  (2)sin  r0,  (4a) 

m,p'  ’  ’  m,p'  '  ^  m,p' 


bP(,)  .  (^1'  lP 

m  Sr^(z)^ 


i  -(1  -  ia(7))r^/r^(2) 


In  Eqs.  (Aa,b),  a  (2)  and  b  (2)  are  complex,  r  (z)  is  the  radiation 

fn  I  p  fn )  p  3 

spot  size,  a(z)  is  related  to  the  curvature  of  the  wavefront  and  is  the 
associated  Laguerre  polynomial. 

Substituting  (3)  into  (1)  and  using  the  orthogonality  properties  of 
L^,  cos  p0,  and  sin  p6,  we  obtain, 

^  y  (  fn,p/  m— l,p  /  m4l,p  y 


fF  (2)1 

.  m,p 

G  (2)  ’ 
V  m,p'  'j 


(5a, b) 


where 


'  (  2  I  '  '1 

A  (z)  =  r  /r  4  i(2m4p4l)  (1  4  a  )c/wr  -  cxr  / r  4  a  /2  ,  (6a) 

m,p  SS  g  ss  ) 

B(z)  =  -  far  /r  4  (1  -  a^)c/«r^  -  a  /2I  -  ifr  /r  -  2ac/wr^'),  (6b) 
\SS  S  y\SS  3/ 


cx'/2], 

(6a) 

2ac/ wr^ j , 

,  (6b) 

*  denotes  the  complex  conjugate  and  the  prime  denotes  a  derivative  with 
respect  to  z,  i.e.,  '  =  3/3z  and 


T  (2)' 
m,  p' 

G  (z) 

.  m,p'  \ 


2n  ® 


(I46  p.)  cos  p0 

p,0' 

sin  p6 


(6c, d) 


4 


where  =  2r^/r^. 

s 

The  function  B(z)  is  arbitrary  and  if  not  specified,  the  equations  for 

a  and  b  in  (5)  are  underdetermined.  If  we  choose  B(z)  =  0,  for 
m,p  m,p 

example,  we  would  in  affect  be  expanding  the  radiation  field  in  the 

conventional  vacuum  Laguerre  modes. ^  Ue  will  show  later  that,  in  general, 

expansion  in  terms  of  the  vacuum  inodes,  B  =  0,  would  require  far  too  many 

modes  to  accurately  describe  the  radiation  beam  over  distances  of  many 

Rayleigh  lengths.  To  find  a  more  appropriate  choice  for  B(z),  we  consider 

the  case  where  the  radiation  beam  at  z  =  0  has  a  Gaussian  radial  profile 

symmetric  about  the  z-axis.  In  this  case  the  complex  radiation  amplitude, 

2  2 

at  z  =  0,  is  given  by  a(r,G,0)  -  q  exp  (-(1  -  ia(0))r  /r^(0))  and  is 

independent  of  0.  Let  us  further  assume  that  for  z  >  0  the  radiation  beam 

profile  remains  approximately  Gaussian  with  a  nearly  circular  cross 

section.  In  this  case  we  expect  the  magnitude  of  the  coefficients,  a  (z) 

m ,  p 

and  b^  ^fz)  to  become  progressively  smaller  as  ra  and  p  take  on  larger 

values.  The  lowest  order  approximation  to  the  radiation  beam  is  then 

given  by  the  a^  q(2)  mode.  From  (5a)  we  find  that  only  the  m  =  0,1  and  p  = 

0  equations  are  relevant  and  are  (0/9z  +  Aq  q)3o  0  ^  ”  ^^0  0  0  =  - 

i(F^  Q  -  Bbq  q).  Ve  now  have  a  specific  expression  for  B(z)  in  terms  of 

one  of  the  moments,  F^  q,  of  the  source  term.  Substituting  B(z)  = 

Fj^  Q(z)/aQ  q(z)  into  (6b)  yields  first  order  coupled  differential  equations 

for  r^  and  a.  Also  using  B(z)  =  F^  Q(z)/aQ  qC^)  allows  us  to  solve  for 

A  in  (6a)  and  hence  for  a  and  b  in  (5a, b). 
m,p  m,p  m,p 
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Ill .  Radiation  Focusing  and  Guiding  in  FELs 


A.  Radiation  Beam  Envelope  Equation 

Ue  first  consider  the  dynamics  of  an  axially  symmetric  radiation 
field  in  the  FEE.  For  a  linearly  polarized  wiggler  field  and  axially 
symmetric  electron  beam  having  a  Gaussian  density  profile,  the  appropriate 
index  of  refraction  for  the  FEE  mechanism^ > ^ 6 » 8 


,  ,  ,  1 
n(r,z,a)  =  1  +  ^  j 


u 


a(r,z) 


(7) 


2  2  2  2  2 

where  w^(r,z)  =  a>j^^(rj^Q/r^(z) )  exp(-r  /r^(z)),  r^^fz)  is  the  electron  beam 

2  1/2 

radius,  r^^  =  ^0  "  ^  I®  r^bO'^^’o^  initial  beam  plasma 

frequency  on  axis,  n^^  is  the  initial  beam  density  on  axis, 


2  . 


a^  =  |e|B^/k^mQC  is  the  normalized  wiggler  amplitude,  B^  is  the  wiggler 
magnetic  field  strength,  k  is  the  wiggler  wave  number,  y  is  the  electron's 
Lorentz  factor,  is  the  electron's  phase  in  the  ponderomot i ve  wave 
potential  and  <  >  denotes  the  ensemble  average  over  all  electrons.  Uith 

the  assumption  that  in  the  source  function  the  complex  radiation  amplitude 
can  be  approximated  by  the  lowest  order  mode,  we  find  that  (2)  can  be 
written  as 


S(£.,z)  =  -^v(a^/r^) 


0,0 


0,0 


(22  'i 


(8) 


2  3 

where  v  =  ~  1^/17x10  is  Budker's  constant  and  is  the 

electron  beam  current  in  :5mperes.  Since  we  are  considering  an  axially 

symmetric  electron  beam  and  radiation  field  we  note  that  a  =  F  =  0 

m,  p  m,  p 

for  p  >  0. 

Using  (8),  equation  (6b)  can  be  used  to  obtain  the  following 
envelope  equation  for  the  radiation  beam, 
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(9) 


r-  .  la^  =  0, 

where 


-4 

(10) 

and  C(z)  =(2w/Y)H(z)a^/ |aQ^Q(z)  1,  H(z)  =(1-F)/(1  +  F)^  and  F(z)  =  r^/rj  is 
the  filling  factor.  The  function  C(z)  measures  the  coupling  between  the 
radiation  and  electron  beam  and  decreases  as  the  radiation  grows-  The 
first  term  on  the  right  hand  side  of  (10)  is  defocusing  and  corresponds  to 
the  usual  diffraction  expansion,  the  second  and  third  terms  are  always 
focusing  while  the  last  term  is  a  defocusing  contribution.  The  envelope 
equation  in  (9)  indicates  that  in  the  high  gain  trapped  particle  regime, 
conditions  for  a  matched  beam  can  not  be  maintained.  However,  in  the  low 
gain  trapped  particle  regime  or  in  the  exponential  gain  regime,  conditions 
for  a  nearly  matched  beam  can  be  achieved.  Using  (5a)  we  find  that  the 
magnitude  of  a^  q(z)  evolves  according  to 

3(rs|ao  q|)/9z  =(4c/w)  ( v/Y)a^r^<sin«i/>/(r^+r^) ,  where  (J^s|ao,ol^^ 

10  2 

proportional  to  the  radiation  power,  P(z)  =  2. 15x10  (  I^q  [Watts] 

B.  Radiation  Guiding  in  the  FEL 

In  the  FEL  the  centroid  of  the  electron  beam  may  be  displaced  off 
axis  by  a  misalignment,  a  redirection  of  the  beam  or  because  of  the 
oscillations  in  the  wiggler  field.  To  determine  the  degree  to  which  the 
radiation  beam  will  follow  or  be  guided  by  the  electron  beam,  we  consider 
the  case  where  the  electron  beam  centroid  is  displaced  transversely  in  the 


=  (2c/w)^ 


-1  4  C^<siniJ/>^  42C<cosiJ/>  4  (a)/2c)r^  C  <sinij/> 


7 


X  direction.  The  index  of  refraction  in  this  case  is  given  by  (7)  with 

w^(r,z)  multiplied  by  (1  +  2(r^X|^/rj^)cos6)  where  is  the  displacement 

of  the  electron  beam's  centroid  and|xj^|<<  r^^.  In  the  FEL  source  term  ve 

consider  only  the  lowest  order  symmetric  and  anti-symmetric  mode  with 

1/2 

respect  to  the  x  axis,  a(r,e,z)  -  (a^  q  1^  cose)exp(-(l-ia)  f,/2) . 

With  this  assumption  the  moments  of  the  source  function,  F  (z),  for 

ni  y  p 

p=0,l,  can  be  evaluated.  For  small  displacements  of  the  electron  beam 
centroid  it  is  easy  to  show  that  the  centroid  of  the  radiation  beam  is 
given  approximately  by 


X  (z)  =  — 1 


where  x,  is  defined  so  that  |a|  is  proportional  to  exp(-((x 

and  (  )n  denotes  the  real  part. 

K 


2,,  2, 
x^)  -  y  )/r^) 


C.  Effect  of  a  Modulated  Electron  Beam 

The  electron  beam  envelope  in  the  FEL  can  undergo  modulations. 

The  modulation  is  symmetric  about  the  z-axis  and  can  be  caused  by  improper 
values  for  the  beam  emittance,  radius  and/or  current  injected  into  the 
wiggler  region. 

For  small  perturbations  about  the  matched  beam  radius,  r^^Q,  we 

find  from  the  electron  beam  envelope  equation  that  r^(z)  =  r^^  (1  +  iysin(KgZ)) 

where  r,  „  =  (2c  /a  k  )  ,  =  a  k  /'I2y  is  the  betatron  wave  number,  due  to 

bO  n  w  w'  ’  B  w  w 

the  weak  focusing  effect  of  wiggler  gradients,  is  the  normalized  emittance, 
and  A  <<  1 . 
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IV.  Numerical  Results 


In  this  section  ve  apply  the  SDE  formulation,  given  by  (5)  together 

with  (6),  to  the  FEL.  Using  the  source  term  given  in  (2)  together  with  (7) 

we  first  present  a  comparison  between;  a)  the  exact  numerical  solution  of 

the  wave  equation  in  (1),  (using  64x6A  Fourier  modes),  b)  the  solution 

using  a  vacuum  Laguerre  modal  expansion  (B=0,  using  10  modes)  and  c)  the 

solution  from  the  Laguerre  SDE  approach  (B  =  F^  0’  modes). 

The  FEL  parameters  used  in  these  illustrations  are  similar  to  those  used  in 

Ref.  9  and  are  given  in  Table  I  where  the  resonant  phase  approximation, 

2 

<exp(-i<i<)>  =  exp(-i«i/j^) ,  is  used  for  demonstration  purposes  and  Zj^  =  nr^(0)/X 
is  the  Rayleigh  length,  X  is  the  wavelength  and  r^(0)  is  the  minimum  spot 
size . 

For  an  axially  symmetric  configuration,  we  show  in  Fig.  2  the  evolution 
of  tha  radiation  beam  amplitude  on-axis  obtained  from  methods  (a),  (b)  and 
(c),  as  a  function  of  propagation  distance.  The  SDE  solution  (c)  is  in 
excellent  agreement  with  solution  (a)  while  solution  (b),  beyond  a  Rayleigh 
length,  grossly  deviates  from  (a)  and  (c).  This  indicates  that  more  modes 
are  required  for  the  vacuum  expansion  solution.  The  excellent  results 
obtained  with  the  SDE  approach  are  also  reflected  in  the  radiation 
amplitude  profile. 

Ve  now  use  the  SDE  method  to  illustrate  guiding  of  the  radiation  beam 

when  the  electron  beam  is  displaced  off-axis.  In  these  numerical 

illustrations,  10  radial  modes  (m  =  0,...,9)  and  2  angular  modes  (p  =  0,1) 

were  used.  In  the  first  example,  the  electron  beam  centroid  is  displaced 

off  axis  according  to  x^  =  x^  ( 1 -sech(k^z ) ) .  Figure  3  shows  the  electron 

and  radiation  beam  centroids,  x.  and  x,  for  x  =  r. /4  =  .075cm  and  X 

b  L  c  b  c 

2ii/k  =  z„/4  =  2.7m.  After  an  initial  transient,  the  radiation  centroid  is 
c  K 


9 


guided  by  and  oscillates  about  the  electron  beam's  centroid.  In  the  next 

example  ve  take  the  electron  beam  centroid  to  be  oscillating  about  the 

2  axis,  Xj^  «=  x^  sin  k^z,  vith  amplitude  x^  =  period  =  Zj^  «= 

10.7m.  Figure  4  shows  the  electron  and  radiation  beam  centroids.  Because 

of  the  high  gain  in  the  radiation  field  the  radiation  centroid  eventually 

follows  the  average  position  of  the  electron  beam's  centroid.  In  the  case 

where  the  electron  beam  centroid  oscillation  is  due  to  the  wiggler  field, 

X  =  a  /vk  and  k  =  2n/X  ,  no  noticeable  change  in  the  evolution  of  the 
c  w  w  c  w 

radiation  field  (compared  to  the  case  for  x^  =  0)  is  observed. 

The  last  illustration  is  for  the  case  where  the  electron  beam  envelope 
is  spatially  modulated.  Using  the  parameters  in  Table  I  ve  find  that  = 
0.06  cm-rad  and  Xg  =  2n/Kg  =  4.66m.  Figure  5  shows  the  amplitude  of  the 
radiation  field  on-axis  as  a  function  of  propagation  distance  when  the 
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Table  I 


Electron  Beam 


Current 

2kA, 

(  V 

=  0.118) 

Energy 

s  = 

50  MeV, 

(y  =  100) 

Radius 

^'bO 

=  0.3 

cm 

Radiation  Bedu. 
Wavelength 
Input  Rover 
Spot  Size 


X  =  10. bum 

A 

P(z=0)  =  BOOMU,  (|a(0,0)l  =  1.84x10  ) 

r^(0)  =  0.6  cm,  =  10.7  m) 


Wiggle r  Field 


Waveleng  t  h 

X  =: 

V 

8  cm 

Wiggler  Strength 

B  = 

V 

2.3  kG,  (a^,  =  1.716) 

Resonant  Phase 

'^R  = 

0.358  rad 
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Figure  Captions 


Fig.  1 

Fig.  2 

Fig.  3 

Fig.  4 

Fig.  5 


Schematic  of  radiation  focusing  and  guiding  in  an  FEL. 


Radiation  amplitude  on  axis,  |a(0,z)l  for  a)  exact  numerical 
solution  (64x64  Fourier  modes),  b)  vacuum  modal  expansion  solution 
(10  modes),  and  c)  SDE  solution  (10  modes)  at  distance  of 
z  =  4z  =  42.8  m. 

K 

Electron  and  radiation  beam  centroids,  x^  and  x^  for  a  displaced 
electron  beam,  x^  =  x^(l-sech(k^z) )  vith  x^  =  r^/4  and 
X  =  2n/k  =  Zj,/4. 

C  C  K 


Electron  and  radiation  beam  centroids,  and  for  an 

oscillating  electron  beam,  Xj^  =  x^sin  k^z  vith  x^  =  t‘j^/4  and 

X  =  2R/k  =  z„. 
c  c  R 


Radiation  amplitude  on  axis,  la(0,z)|  for 
beam,  r^  =  +  6  sin(KgZ))  vith  = 
Xg  =  2n/Kg  =  4.66  m. 


a  modulated  electron 
0. 3cm,  6=0.1  and 
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Appendix  V: 

Radiation  focusing,  guiding  aiul  steering  in  the  Free 
Electron;  Laser 


1:5 


Radiation  foclsing.  glmding  and  steering  in  free  electron  lasers 

P.  Spt»n|l»,  A.  Ting*.  B.  Hnliil**  and  C.  M.  Tang 


Naval  Krsrarch  Laboratory 
Vashinglon,  DC  20375-5000 


Abs  tract 

In  a  !r»»  »l»c:ron  laser  (FED.  'he  radiation 
(•.eld.  wiggler  (leld  and  electron  beaa  resonantly 
couple  and  aodify  the  refractive  indev  in  the  vicinity 
of  the  electron  beaa.  The  refractive  inde*  Is 
modified  such  that  the  radiation  beam  will  tend  to 
foc'js  upon  the  electron  bean.  A  method  for  solving 
the  3-D  wave  eouation  for  the  FEL  process  is  outlined. 
This  approach,  called  the  source  dependent  expansion 
method,  provides  an  excellent  analytical  and  numerical 
technique  lor  studying  optical  focusing,  guiding  and 
steering  in  FELs.  A  radiation  envelope  equation  is 
derived.  Conditions  and  parameters  necessary  to 
achieve  guided  radiation  beans  (constant  radius)  in 
the  exponential  gain  regime  are  obtained  for  FELs 
driven  by  either  induction  linacs  or  rf  linacs. 
Imieedlately  prior  to  saturation  in  the  exponential 
gain  region.  the  ponderomotlve  potential  is  large 
enough  -o  trap  the  beam  electrons.  The  viggler  field, 
at  this  point.  could  be  tapered  to  further  increase 
the  operating  efficiency.  The  possibility  of  bending 
or  steering  radiation  beams  in  FELs  is  discussed  and  a 
condition  necessary  for  radiation  guiding  along  a 
curved  electron  bean  orbit  is  obtained. 

Introduction 


In  many  short  wavelength  free  electron  laser 
devices  the  radiation  bean  will  not  be  confined  or 
guided  by  a  structure  such  as  a  waveguide. 
Fur  t  hernor  e .  in  order  to  provide  high  gain  and 
efficiency.  it  is  usually  necessary  for  the 

interaction  length  (length  of  wiggler  field)  to  be 
long  compared  to  the  diffraction  length  (Rayleigh 
length)  associated  with  the  radiation  beam.  In  the 
FEL  the  tendency  of  the  radiation  beaa  to  diffract 
away  over  a  distance  of  a  few  Rayleigh  lengths  can  be 
overcome  by  a  focusing  phenomenon  arising  from  the 
resonant  coupling  of  the  radiation  and  viggler  fields 
with  the  electron  beam  (1.21.  This  radiation  focusitig 
effect  plays  a  central  role  in  the  practical 
utilization  of  the  FEL.  This  phenomenon  was  first 
analvzed  for  the  low  gain  FEL  with  transverse  effects 
where  it  was  shown  that  the  diffractive  spreading  of 
the  radiation  beam  could  be  overcome  by  a  focusing 
effect  arising  from  the  modified  index  of  refraction 
I ’.  1  .  Optical  guiding  in  FELs  operating  in  the  small 
signal  exponential  gam  regime  has  been  studied  for 
the  asvmptotic  behavior  of  the  radiation  beam  |3-6). 
Recent  Iv.  a  general  formalism  for  optical  focusing, 
guiding  and  steering  has  been  developed  and  applied  to 
FELs  I  7), 

In  the  following,  we  employ  a  modal  expansion 
technique  to  examine  the  optical  beam  as  it  propagates 
'hroitgh  the  wiggler.  The  formalism  has  the  merit  that 
with  only  a  few  modes  it  permits  an  accurate  solution 
of  the  wave  equation  throughout  the  interaction 
region. 


Model 

In  our  model,  'he  vector  potential  of  an  axially 
svmmetric.  linearly  polarized,  radiation  field  is 

Ap(t.z.t)  .  A(r.z)  e‘^““'''“' U^/2  •  c.c.,  (1) 

^(r.z)  Is  fh#*  rowplpx  radiation  field  amplitude, 
u)  is  rhe  frequency  and  c.c.  denotes  the  coaplex 

on  HI jfa  f  e  , 


rhe 


va  ve 


equation  governing; 


i  s 


I  ’-fr 

L.] 

r  3r  i 

3r  j 

(2) 


where  J  (r.z.t)  is  the  driving  cuiient  density 
Substituting  (1)  into  (2)  leads  'o  the  following 
reduced  wave  equation. 
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a{ r .z ) 


S( r . z .a  ) . 


(3) 


vhere  a(r,2)  •  -  |a|exp(i4)  is  the 
noraalized  coaplex  radiation  aaplitude  and  we  have 
assuaed  that  a(r.z)  is  a  slovly  varying  function  of  7. 
i.e.,  |(}a/3z)/a|  <<  im/c.  The  source  function.  S.  is 
given  by. 


2ii/u 

-  ^  r  j 

c  J  > 
o 


( r . z . t  )e 


i  (luz/c  -<vt  ) 


d  t  . 


(4  ) 


It  is  possible  to  relate  the  source  function,  S. 
to  the  index  of  refraction  associated  with  the  medium 
by  noting  that  the  wave  equation  for  A^  in  a 
nonmagnetic.  t ime- independent ,  nonlinear  medium  is 

(  -  (n^(  r  .  z  .a  >70^ )  /  3t  ^  )i^  =  0.  where  n  is  the 

index  of  refraction  associated  with  the  medium  and  is. 
in  general,  complex  and  a  nonlinear  function  of 
a(r.z).  Comparing  the  reduced  wave  equation  wtitten 
in  terms  of  n(r.z.a)  with  (3)  we  find  that  the  source 
function  can  be  written  in  terms  of  n. 


S(r,z,a) 


((>i/c)^|l-n^(r,z.a)ja(r.z). 


(5) 


Source  Dependent  Expansion  Method 


In  order  to  solve  (3)  we  will  use  the  source 
dependent  expansion  (SDE)  method  |7|.  In  this  method, 
we  choose  the  following  representation  for  a(r.z)  in 
terms  of  Laguerre-Gaussian  functions. 


a(r.z)  .  Y.  *B^**''m 


lr‘(z) 


-(l-ia(z))r  'r  (z) 
e  ^  .  (6) 


vherv  ■  •  0»1.2.‘**.  In  Eq .  (6),  a  (2)  are  the 

complex  amplitude  coefficients,  r  (z)  is  tRe  radiation 
spot  size.  «(z)  is  related  to  the  radius  of  cuAvature 
of  the  radiation  beam  wavefront.  R  -  -  (uw2c)D'a 

and  L  is  the  Laguerre  polynomial.  Solving  for  the 
unknown  quantities  a  ,  r  and  a  in  terms  of  the  source 
term  S  allows  us  to'com^letely  describe  the  radiation 
dynamics.  The  representation  in  (6)  is  underspecified, 
since,  when  (6)  is  substituted  into  (3)  and  moment>:  nf 
the  source  function  taken,  there  remain  moie  unknown 
quantities  than  available  equations.  The  additional 
degrees  of  freedom  in  our  representation  allow  us  tn 
specify  a  particular  functional  relationship  lot  i  hv 
unknown  quantities  r  and  a  in  such  a  wav  that,  if  ilir 
radiation  beam  profiJe  remains  approximately  baussiau, 
the  number  of  modes  needed  to  accurately  desciibe  the 
radiation  beam  is  small.  This  yields  the  following 
first  order  coupled  differential  equations  for  and 


r^  -  2co/ur^  .  -  r^  . 


2(1' 


^ 

a  )  c  /  ur 

s 


2(H„ 


oHj  !  . 


(  ?h  I 


•nH  j  of  first  order  ordinary  dift»renel«l 

eouATions  lot  the  coeplen  jepllludes  *^(1), 

*-  •  •  -  'K  ■  <^<=> 

where  fl  •  ^1  *  '  "  3/Jt,  *nd  (  )  denotes  the  reel 
end  laa^inAfy  ^ert  of  the  enc losed * lunc 1 1  on .  In  Eqs- 
f^).  the  functions  A^.  B.  end  ere  Jiven  by 

»  r’'r  *  -  or'^'r  •  o'/*]. 

'  s  S  S  S  S  / 


^  ( z)  •  lelA  /  rm  c(e*p(lk  r)  •  c.c.l  e  T. 
-V  '  V  o  V  » 


Is  the  wlgf^le  velocity,  r  is  the  Lotemr  fertor.  A  j 
the  vector  potentlel  eoplltude  of  the  vicgle* 

field  end  k  •2t/x  Is  the  wiggler  vAve  ou^hor 

Substituting  ^6)  inlfo  the  expression  (or  S.  (4),  gives 


S- 


JT) 


7iL/m 


-  i 


t  -  T(  z  .  t  ) 
O 


or  ■  r  •  f  1  -  o‘  )c 
I  s  s 


/r  - 
s  s 


2«c/ 


F^(2)  ■  i  S(  C,  2  )L^(  Oexpt  -  (  1  *  ia)  C'T  ) , 
o 

where  C  -  2r^/r^. 

s 

The  aerits  of  the  SOC  aethod  can  be  deaonstrated 
in  a  coaperison  between;  a)  the  exact  nuaerical 
solution  of  the  wave  equation  in  (3).  (using  64x6^ 
Fourier  aodes).  b)  the  solution  using  a  vacuua 
Uaguerre-Gaussian  aodal  expansion  (10  aodes)  and  c) 
the  solution  froa  the  La^erre-Caussian  SDE  approach 
(10  aodes).  Figure  1  shows  the  radiation  beaa 
aaplltude  on-axls  obtained  froa  aethods  (a),  (b)  and 
(c)  after  four  Rayleigh  lengths  for  the  FEL  paraaecer 
in  Table  I-  The  SDE  solution  (c)  Is  in  excellent 
agreeaenr  with  solution  (a)  while  solution  (b),  beyond 
a  Rayleigh  length,  grossly  deviates  froa  (a)  and  (c). 


Fig.  I  Radiation  amplitude  profile.  |a(r.z)|  for: 

a)  exact  nuaerical  solution  (6^x64)  Fourier 
aodes).  b)  vacuua  aodal  expansion  solution 
(10  aodes).  and  c)  SDE  solution  (10  aodes)  at 
a  distance  of  z  • 

Refractive  Index  Associated  with  FELs 

In  the  following  derivation  of  the  refractive 
index  associated  with  the  FTL,  a  nuaber  of  siapllfylng 
assumptions  are  aade.  Ve  assuae.  for  exaaple.  that 
the  beaa  electrons  are  aonoenerge t ic  without  betatron 
oscillations  and  that  the  radiation  is  of  a  single 
frequency  |8|.  To  obtain  an  expression  for  the 

refractive  index  we  write  the  nonlinear  driving 
current  density,  J  .  as 

^  -X 

J  •  -|e|n.(i)v  (z)v  fi(2-r(t.t  ))dt..  (81 

-X  '  b  -w  02j  o  0 

where  n.  fr)  is  the  ambient  beam  density,  v  ts  the 
axial  electron  velocit^  at  z  ■  0.  is  tnf  tiae  a 

given  electron  '•rosses  the  z  ■  0  plane. 


where  a  •  lelA/ac*'.  T*  t  •  [dz'^v(r'.:  )anf) 

w  w  o  o  J  z  o 

o 

the  t  integration  is  over  all  entry  times.  Equating 
(^)  w?ih  (5)  and  carrying  out  the  integration  over  t. 
we  find  the  index  of  refraction  associated  with  the 
FEL  to  be  given  by 

nj^j(r,2,a)  -  1  ♦  (M^(r)/2«^)  ^  .  (10) 

where 

z 

♦  -  I  '  iln(a/|a|)  -  v^(  z .  )  ]dz  * 

0 

Is  the  relative  phase  between  the  election  and  the 
ponderoao t i ve  wave,  initial  phase  of 

2n 

a  given  electron  and  ^  “  C2ii)  ^  J  is  an 

^  '  0 

0 

^nseablc  average  over  the  initial  phases.  The  radia; 
profile  of  the  index  of  refraction  as  given  by  to. 
(10)  supports  self- focusing  of  the  radiation  in  an 
FEL.  It  should  be  noted,  for  conpleteness.  that  the 
relative  phase  satisfies  the  penduiui*  equation  given 
by 

-  t'^(i.i/c)  [4Ja^/3z  sin.j.  (11) 


Radiation  Bean  Envelope  Equation 

Equations  (7a)  and  (7b)  can  be  combined  -o  give 
the  following  envelope  equation  for  the  radiation  beam 


K  r 


0. 


(12) 


where 

K^-(2c/»)‘ (-l‘C^<sin»>‘'.2C<cos*>.(u/2c)r‘C  Csin*.'|r^  . 

(U) 

and  C(z)  •  (2v/T)G(z)a  /  |a  (z)l,  aeasutes  the  coupling 
between  the  radiatioR  and  electron  beam, 

-1  1 

w  .  r|^'2c)‘  -  Ij^/17xl0'  is  Btidker's  constant, 

I.  Is  the  electron  beam  current  in  amperes, 

U  J  1 

G(z)  .  (l-f)/(l.f)  and  f(z)  .  (r./r  )‘  is  the  filling 

b  s 

factor  associated  with  a  Gaussian  electron  beam 
density  profile.  The  first  term  on  the  right-hann 
side  of  (13)  is  the  usual  diffraction  teim,  the  second 
and  third  terms  are  focusing  while  the  las!  ’ei- 
provides  a  focusing  or  defornsing  cont r i btn  i on .  In 
the  high  gain  trapped  particle  legime.  <sin*>  and 
<cose>  are  approximately  constant,  while  |a  (:)l 
Increases  with  z.  Hence.  K  depends  on  z  and  a  guided 
beam  (r'  •  0)  cannot  be  exactly  maintained  in  ihis 

regime.  although,  the  radiation  envelope  is  still 
reasonably  we  1 1  -  con f i ned  .  In  the  low  gain  iiapped 
particle  regime  {a  (z)|  increases  slightlv  am 
therefore.  a  guide?  beam  can  be  appr o» i ma t e  1  v 


In  e»t^'er  rhe  Cc^pton  oi  ftjvAn  *«Ponentl«l 
f»in  conditions  for  «  stabi»  fuid«d  b*%m  can 
be  found. 


Cuidfd  Radiation  Bfaas  in  ihf  Enponantial  Cain  Kfulag 


In  this  section.  we  obtain  the  necessary 
conditions  to  achiew*  yuided  radiation  beaas  in  both 
the  Coepton  (noncollective)  and  Kasan  (collective) 
exponential  gain  regiaes.  Bv  considering  the  lowest 
order  aode  (Canssian  profile)  we  find  "hat  the  source 
ter«  appropriate  (or  the  high  gam  Coepton  and  Raaan 
regiee  is.  respectively, 


1 


(wi,(r)/c)'(a  it  f,)* 

S(  r .  z  )  • - c — -  a<  r  .  r  ) 

Y(l*a;.'2) 


(ak-iD' 

1/2 


(r  Xfik-iD 


(Ua.b) 


where  6k  and  F  are  the  wave  nuaber  shift  and  growth 
cate  respectively  and  f.  is  the  usual  difference  of 
Bessel  functions  due  to  the  linear  wiggler.  The 
lowest  order  aode  Is  taken  to  have  the  fora 

a(r.i)  -  a^(0)exp(ij(ak-irid2' -d-ialr'/rp.  (15) 

o 


For  rhe  purposes  of  : 1 lus t ra t ion .  ve  will  consider  the 
Coepton  FEL  regime  in  which  the  electron  bea«  has  a 

•1  ^ 

Gaussian  densitv  profile.  n.(r)  •  n  exp(-  r*  rT). 

DO  O 

The  conditions  for  a  guided  radiation  beaa  require 
that  the  waist  and  curvature  of  the  radiation  beaa 
reaain  constant.  <r'  .  a'  .  0).  Setting  r'  «  a' .0  in 

Eqs.  {2a. b)  and  solving  for  T,  6k.  r  ,*and  a.  the 
folloving  results  for  a  guided  beaa  are  obtained. 


r  .  (i,.a‘)‘*{!.2f)''r 


6k. aT. 


(16a. bj 


fll 

'  V/ 


1/4 


,2/4  ,  I, 

2  srfj 


1/2 


(K3f/2) 


.(16c) 


r  (f-1)  .  0.25  X 

S  V  \\tj 


(l-a;/2)^''" 

,  i/2  1/2' 


(16d) 


(f'  (:.3f )) 


1/2 


(16e) 


1/2  .1/2  -t  -* 

where  =2f_(u/v)  a  V  (l.a'’/2)  and  fvr^'r'  is 

0  B  WWW  b  s 

the  filling  factor. 


Figure  2  shows  the  spatial  evolution  of  the 
radiation  waist  for  the  Induction  linac  driven  FEL 
paraaeters  In  Table  I.  The  pacaaeters  in  Table  I  are 
consistent  with  Eqs.  (16)  and  have  been  chosen  to 
produce  a  guided  radiation  beaa  in  the  Coapton 
exponential  gain  regiae.  The  guided  bea*  conditions 
can  be  shown  to  be  stable  |9|,  this  is  shown 
nua^rically  by  changing  the  spot  size  of  the  injected 
radiation  beaa.  Figure  3  shows  that  irrespective  o( 
the  initial  value.  the  spot  size  asyaptotes  to  the 
aatched  (guided)  beaa  value.  Figure  4  shovs  the 
evolution  of  the  spot  size  for  the  rf  llnac-drlven  FEL 
paraaeters  in  Table  II.  As  in  Table  I.  the  paraaeters 
in  Table  11  have  been  chosen  to  produce  a  guided 
radiation  heaa  tn  the  Coapton  exponential  gam  regiae 
and  are  consistent  with  Eqs.  (16). 


Fig.  2  Spatial  evolution  of  the  radiation  spot  size 
n  the  exponential  gain  regiae  for  Induction 
linac  driven  FEL  paraaeters  given  in  Table  I. 


Fig.  3  Spatial  evolution  of  the  radiation  spot  size 
in  the  exponential  gam  regiae  for  initial 
spot  sjzes;  a)  0.35  c»,  b)  0.24  cn.  and  c) 

0. 15  c». 


Fig.  4  Spatial  evolution  of  the  ladiation  spot  size 
In  the  exponential  gam  regme  (ot  rf  linac 
driven  FEL  parameters  given  m  Table  TI 


Frrr  tl»ctron  driven  by  either  induction 

or  rf  linic?  could  initially  operate  in  the  fuided. 
evponential  gain  regiae  until  saturation  occurs. 

prior  to  saturation.  the  ponderoaot  Ive 
potential  can  be  large  enough.  as  in  the  above 
illustrations,  to  trap  a  signillcant  fraction  of  the 
electrons.  At  this  point,  the  uiggler  field  can 
be  spatiallv  tapered  to  achieve  a  significant  Increase 
In  the  operating  efficiency  and  a  soaewhal  saaller 
input  signal  into  the  FEL  aap] i f ler  . 


To  deteraine  the  viability  of  tapering  the 
viggier.  prior  to  saturation,  the  trapping  potential 
associated  vith  the  ponderoaot ive  uave  is  needed.  For 
llnearlv  polarized  waves.  the  fractional  trapping 
potent lal  is 


2 

r«  c 
o 


2/2 


(17) 


The  radiation  aaplicude  at  saturation  can  be  obtained 
(roe  the  intrinsic  efficiency  of  the  FEL.  Using 
arguaents  based  on  electron  trapping  in  the 
ponderoao  t  i  ve  wave,  we  find  that  the  intrinsic 
efficiency  in  the  exponential  (aaxiaua)  gain  regiae  is 


n 


(18) 


Using  the  induction  llnac  paraaeters  in  Table  I  as  an 
illustration,  ve  find  that  the  intrinsic  efficiency  is 
h  •  •  0.66Z.  Froa  this,  the  fractional 

trapping  potential  at  the  end  of  the  exponential  gain 

regiae  is  je | a^^^^/va^c^-bS.  aaxing  it  possible  to 

trap  the  electron  beaa  while  tapering  the  wiggler 
field.  In  addition.  the  Initial  fractional  energy 
spread  of  the  electron  beaa  aust  be  soaewhat  less  than 
n.  This  places  a  liaitatlon  on  the  fractional  energy 
spread  of  the  electron  beaa,  ^  0.66Z.  One 

contribution  to  the  beaa  energy  spread  is  the 

transverse  eaittance,  &E/E.  .  (l/2)(t  /r,)^. 

0  no 

Therefore,  the  noraalized  bee i  eaittance  aust  satisfy, 

<  {2(*/k.^)'''^r.  .  0.034  ca-rad. 
n  ^  D 


and  Guiding  of  Radiation  Beaas 


Using  the  SOE  foraallsa,  it  is  possible  to 
discuss  the  bending  of  a  radiation  beaa  by  a  curved 
electron  beaa  in  an  FEL.  For  saall  displaceaents  of 
the  electron  beaa  centroid,  a  nonaxisyaaef ric  aodal 
expansion  slailar  to  (6)  can  be  perforaed  and  the 
spatial  evolution  of  the  centroid  of  the  radiation 
beaa  found.  Figure  5  shows  the  centroids  of  the 
electron  and  radiation  beaas  for  an  FEL  in  the  trapped 
particle  regime  with  paraaeters  given  in  Table  I. 
Steering  of  the  radiation  beaa  by  t’  electron  beaa  is 
clearly  deaonslialed  in  this  flguri 


It  Is  interesting  to  consider  the  conditions 
under  which  the  radiation  beaa  could  be  guided  by  a 
curved  electron  beaa,  as  shown  in  Fig.  6.  Such  a 
situation  could  aake  possible  a  cyclic  FEL  driven  by, 
for  exaaple.  .a  betatron  generated  electron  beaa.  In  a 
cyclic  FEL,  the  radiation  beaa  would  be  guided  by  a 
circular  electron  beaa.  The  wiggler  field,  which  is 
along  the  circular  orbit  of  the  electron  beaa,  cannot 
be  spatially  contoured.  Therefore,  in  the  trapped 
particle  tegiae,  enhanceaent  of  the  FEL  efficiency 
aust  be  achieved  bv  inducing  an  accelerating  electric 
field  along  the  beaa  orbit.  For  cyclic  electron 
beaas,  the  induced  electric  field  can  be  generated  by 
increasing  the  aagnetlc  flux  within  the  orbit  of  the 
electro."!  beaa. 


Fig.  S  Electron  and  radiation  beam  centroids,  x 
and  Xj^  for  a  displaced  electron  beam^ 
"h  ■  x„(l-sech(k  z))  with  x  .  r./c  and 
X“  .  2ii/6  -  4Z„. 

c  c  R 


Fig.  6  Configuration  showing  guiding  of  radiation 
hesm  by  a  curved  electron  beam  with  radius  of 
curvature,  R  . 

0 


To  exaaine  the  conditions  under  which  guiding  can 
be  achieved  in  the  exponential  gam  regiae,  ve  denote 
the  radial  position  by  r  =  R  »  x.  where  R  is  the. 
radius  of  curvature  of  the  election  beaa  and  x  is  the 
radial  displacement  from  the  center  of  the  curved 
electron  beaa  (see  Fig.  6).  The  FEL  refractive  indev 

(correct  to  order  x/R  )  is 
o 

n  •  n,  ,  ♦  x/R  ,  ( I'll 

fel  0  1  ‘  ' 

where  n.^  is  given  by  (10).  In  the  exponential  gain 
regiae.  S‘guided  radiation  beam  in  a  curved  FEL  is 

possible  if  R  >  R  ,  where 
0  -  Bin 

"ain  ■  ^'l''"*''’fel^l-  <■'” 

Substituting  the  expressions  foi  F.  aiui  i  .  f  i  i,m 
Eqs.  (16),  into  (20)  yields 


4(  1  ,  Ofv'  I. 

g  _  _ 2 _ _ 

(1.2f)(3f.2)^^'fga^(l.a,;  '2  ) '  ' '' (  v/ , ) " '  ‘  ' 


R  ,  (f-1) 

■  in 


'Tib' 


f  ot  a  rin»e  r  :  c  a  i 


consider  the 


''  T  j  on 


e«a«plr  of  ^  consider  the 

■  in 

following  pataaeters.  y  •  100.  I  •  T  k> .  r.  •  0.3  c». 

o 

•  ^  •  I  . 'r .  !  !  *nd  .  0.0^  fT»bl»  !>  Fot  th*s» 

p*r»»»i»rs.  'hf  lurning  radius  required  for  « 

guidrq  radiation  t>«a>  is  R  •  4^^  a. 


Table  I 

Paraae’ers  Associated  with  an  Inauction  linac  3ri«en 
FEL  in  ;ne  Exponential  Tain  Regiae 


Electron  Beaa 
Current 
Ener  g-y 
Radius 
Eai t  tance 

Vigglecfield 
Wavelength 
Vigglec  Strength 


I  .  21t>.  (  \)  .  0.  1 18) 
E^  .  50  MeV.  ( T  .  100) 
-  0. 3  CB^ 

c°  <  34x10’  cB-rad 
n 


X  .  8  CB 


2.3  ItC  (a 


Radiation  SeaB 
Wave  1  eng  t  h 
Spot  Sue 
'guided  beaa) 
e-£olding  length 


X  .  10.6  UB 

r  .  0.:i  CB.  (2 
s 

L  .  1/r  .  64  ca 


Intrinsic  Efficiencv 


Saturated  Pover 


Trapping  Potential 


P  .  660  hU  (a  .  7x10  ) 

sa  t 


e  ♦. 

trap 


c‘  .  6.0J 


Table  II 

Paraaeters  Associated  wi th  an  RF  Linac  Driven 
FEL  in  the  Exponen tial  Cain  Regiae 


Electron  Beaa 
Peak  Current 
Energy 
Radius 
Eai t  tance 

Higgler  Field  (planar) 
Wavelength 
Higgler  Strength 


I.  .  500  A 
E°  .  150  MeV 

0  a 

ty,  •  1  , 

c  <  7x10'  CB-rad 
n  - 


X  -  ::  cm 
8*'  .  ‘*00  C  (a 


Radiation  Beaa 
Wavelength 
Spot  Size 
(guided  beaa) 
e-folding  length 


r^(0)  .1.1 


Intrinsic  Efficiency 


n  .  4k/k  =  0.25Z 


Saturated  Power 


^sat  ’  KH  (a  .  7.25x10  ) 

l'l^rap'^o'^  ■ 


Free  electron  lasers  di  i  ven  hv  .|>hrr  mdu 
llnacs.  such  as  the  ATA,  nr  high  cuiT.ni  r*  iinx. 
operate  In  the  guided,  exponential  gain  legi*.  -jittil 
saturation  occurs.  At  this  point,  i  he  viggjer  fieici 
could  be  spatially  tapered  so  a»  to  opetxte  -he  FEE  in 
the  trapped  particle  regime  in  order  rn  fnrthei 
increase  the  operating  efficiency 

He  also  exaained  the  possihiliiv  of  bending  ni 
steering  radiation  bea»s  in  FELs  We  find  a  -ondi'inn 
which  places  a  lower  lieit  on  'he  'adt"S  of  -'it.x'-.nn 
of  the  electron  bea*  necessary  for  the  radix' ion  :o  be 
guided  along  a  curved  path. 
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Conclusion 

The  source  dependent  expansion  (SDE)  Bethod 
provider  an  excellent  analytical  and  nuBerlcal 
technique  for  studying  optical  focusing,  guiding  and 
steering  in  "ELs.  He  find  that  guided  radiation  beans 
in  the  FEL  can  be  achieved  both  In  the  CoBpton  and 
Reman  exponential  gam  regimes  but  cannot  be 
Baintained  in  the  high  gain  trapped  particle  (tapered 
wi gg ler  )  regime. 


Appendix  VI; 

Optical  gain,  pliase  shift,  and  profile  in  tlie 
Free  Electron  Laser 


\\ 


I.  Introduction 

A  vell-kjiovn  feature  of  the  f ree-elec t ron  laser  (FEL)  is  that  the 

refractive  index  of  the  nediun:  is  a  complex  function  and  hence  the 

radiation  is  amplified  and  to  some  extent  focused  in  the  vicinity  of  the 

1  2 

electron  beam.  ’  It  may  then  be  possible  for  the  electron  and  radiation 
beams  to  interact  over  an  extended  length  along  the  viggler,  with  the 
diffractive  tendency  being  compensated  by  the  PEL  interaction,  thereby 
enhancing  the  efficiency  of  the  process. 

Considerable  progress  has  been  made  in  studying  this  process  by 

3-8 

several  authors.  The  purpose  of  this  paper  is  to  apply  the  formalism  of 
the  Gaussian-Laguerre  modal  source  dependent  expansion  (SDE)  of  Ref.  8  to 
examine  the  propagation  and  guiding  of  the  optical  wave  in  an  amplifier 
operating  in  the  exponential  gain  regime,  for  a  variety  of  operating 
condi tions . 

The  plan  of  this  paper  is  as  follows.  In  Section  II  the  formalism  of 
the  SDE  is  employed  to  obtain  the  evolution  equations  for  the  radius  and 
the  curvature  for  the  lowest  order  mode  of  the  optical  beam,  along  with  the 
relevant  dispersion  relation  for  a  Gaussian  electron  beam  driving  an  FEL 
amplifier  in  the  small  signal  regime.  In  Section  III  numerical  solutions 
of  the  single-mode  equation  for  the  radius  of  the  optical  beam  are 
presented  and  compared  to  the  result  from  a  multi-mode  truncation  of  the 
radiation  field.  In  this  case,  and  for  cases  not  presented  herein,  the 
single-mode  and  multi-mode  results  indicate  that  the  radiation  beam  profile 
entering  the  viggler  asymptotes  to  a  unique  form  after  an  initial 
transient.  Additionally,  the  numerical  values  of  the  radius  of  the 
radiation  envelope  and  of  the  wavefront  curvature  are  in  fair  agreement, 
irrespective  of  the  degree  of  mode  truncation,  indicating  the  usefulness  of 
the  single-mode  equations.  Limiting  ourselves  to  these  equations,  the 
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electron  beam  Is  then  allowed  to  oscillate  at  the  betatron  wavelength  and 
the  resulting  radiation  profile  examined.  It  is  found  that  the  optical 
beam  envelope  follows  that  of  the  electrons  with  almost  identical 
wavelength,  but  retarded  in  phase.  Section  IV  discusses  the  results, 
deriving  formulae  for  the  matched  radiation  beam  profile  (i.e.,  radius  and 
curvature)  in  terms  of  the  electron  beam  and  viggler  parameters.  It  is 
shown  analytically  that  perturbations  of  the  profile  are  spatially  damped 
out,  consistent  with  the  numerical  observations  indicating  a  unique, 
asymptotic  matched  radius  and  curvature.  Appendix  A  presents  the  necessary 
details  required  to  derive  the  source  term,  for  the  wave  equation,  for  a 
planar  viggler  and  an  electron  beam  with  uniform  density  along  the 
direction  of  propagation.  Appendix  B  considers  the  effect  of  the 
modulation  of  the  electron  beam  on  the  optical  wave.  Specifically,  a 
simple  analysis,  taking  into  account  sideband  generation,  indicates  that 
the  dispersion  characteristics  of  the  primary  wave  are  only  slightly 
modified  for  typical  experimental  parameters.  Appendix  C  presents  the 
details  of  the  stability  calculation. 


II.  Mathematical  Formulation 

The  purpose  of  the  present  section  is  to  present  the  salient  features 

g 

of  the  source  dependent  expansion  method  so  as  to  fix  the  notation  and  for 
reference  in  the  subsequent  sections. 

For  a  planar  viggler,  it  is  appropriate  to  assume  a  linearly  polarized 
radiation  vector  potential 


A  =  (  1/2)  A(r,e,z)exp 


■ 

'OJ  z 

i 

-  -  <0  t 

c 

. 

. 

e  *  c.c. , 

-X 


2 


vilh  angular  frequency  w  and  complex  amplitude  A.  In  the  slovly  varying 
envelope  approximation,  the  vave  equation  reduces  to 


13  3  13^  Oa  S(r  6 

r  d  D  j 


where  a  -  |e|A/n^c^,  and  the  source  function  Is  given  by 
s(r,e,z)  -  -  (j  (r.e.Oexp  -1  —  -  w  t  ] 

B  I  LI 


(1) 


(2) 


Here  e  is  the  charge  on  an  electron  of  (rest)  mass  m  ,  J  (r,0,r)  is  the 
current  density  and  {  ^slov  that  only  the  spatially  and 

temporally  slov  part  of  the  quantity  in  braces  is  to  be  retained. 

The  basic  premise  of  the  vork  presented  herein  is  that  the  radiation 
field  is  azimuthally  symmetric  and  the  vector  potential  is  expressible  as: 


a  (r,e,z)  =  Yj  (3) 

m=o 

vith  D  =  L  (  £,)  exp  ( -  [  1  -  i  a(  2 )  11/2 )  ,  where  ?:  =  2r^/r^(z),  r  (z)  is  related 

SD  ID  S  S 

to  the  radiation  spot  size,  a( z )  is  proportional  to  the  curvature  of  the 

wavefront,  and  L  ( ^)  is  the  Laguerre  polynomial  of  order  m. 
m 

Nov,  if  the  transverse  profile  of  the  radiation  beam  is  close  to  a 

3  5  7  8 

Gaussian,  the  lowest  order  mode  is  expected  to  dominate  ’  ’  ’  ,  and, 
following  Ref.  8,  it  is  simple  to  show  that  the  associated  vector  potential 
evolves  according  to 


and  the  spot  size  ind  wavefront  curvature  evolve  via 
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-  nu-"j  -  “U-JiJ 

u  r  OK  o 


vhere 


A  -  d  •  f/i  2,  c  ad  Id] 

o  r  dz  s  L  2  r  dz  s  2  dz  J  ’ 


(or  s 
s 


the  F's  are  given  by  the  following  overlap  integral: 


-I;;!  « s(5,.)D„‘(t,z), 


and  the  label  R  (I)  indicates  the  real  (imaginary)  part. 

Noting  that  =  1*  the  normalized  vector  potential  is  seen  to'be 

given  by  (Eq.  (3)] 


a(r,e,2)  =  a  (z)  exp<^-(  l-ia(z) 

I  r^(z)J 

vhere,  in  the  exponential  gain,  small-signal  regime, 


^(z)  =  a(o)  exp  |i  dz^  Ak(z^)-ir(z^)  i. 


Bere  a(o)  is  the  input  signal  at  z  =  o,  and  the  two  components  of  the 
refractive  index  are  given  by 


, .  c  Ak  .  .  c  2  S  1  -  i  o'! 


Assuming  the  electron  beam  profile  to  be  given  by 


-  n 


bo 


’’bo 

^  1 

r  2 

r 

rb(r) 

exp 

r^(2). 

(10) 


where  is  the  electron  beam  radius  at  r  and  n^^^  is  the  beam  density  at 

r^(z)  -  r^^,  the  source  term  in  Eq.  (1)  may  be  readily  evaluated  (Appendix 
A),  to  obtain 


S(r,z) 
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‘^o 


■B  .  3  2 
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(  2  ) 
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i:b(z) 

exp 

2 

^  "b 

u  Ic  a  a 

V  V 

c(dk-ir) 
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where  the  vector  potential  of  the  planar  viggler  of  periodicity  2ii/k^  is 
given  by 


A  =  A  cos(k  2)e  , 

-V  V  v  -x 


(12) 


a  =  e  A  /m  c  , 

V  V  0 


(13) 


T  is  the  relativistic  mass  factor,  f-  is  the  usual  difference  of  Bessel 

D 

functions,  f  =  J  (Q  -  J.(C),  C  «  (l/4)aj/[l  *  (l/2)aj],  and 

D  O  1  V  V 


“bo  - 


1/2 


is  the  plasma  frequency  of  the  electron  beam  with  density  . 

Substituting  Eqs.  (8)  and  (11)  into  Eq.  (6)  and  making  use  of  Eqs.  (M 
and  (5),  it  is  simple  to  show  that  the  equations  reduce  to 


d  a _ 

d(k^z) 
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(k  r  )‘ 
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(14b) 
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0,  (Uc) 


vhere 


B  Kck 


r^(z) 


V  l  b  s)  .  I^y 

2y^  r  21^ 

"  k2(r,/rp2] 


(Ud) 


The  spatial  evolution  of  the  system  is  governed  by  the  differential 
system  (14a)  and  (14b)  along  vith  the  dispersion  relation  (14c),  the 
solution  of  vhich  yields  a(2),  r^(z),  tik(z)  and  r(2). 


III.  Numerical  Results 

Having  obtained  the  single  mode  system  of  Eqs.  (14),  it  is  of  interest 
to  determine  the  extent  to  vhich  it  approximates  the  general  solution  in 
(3).  Once  it  is  established  that  Eqs.  (14)  provide  an  adequate 
representation  of  the  general  solution,  it  is  then  possible  to  study  a 
variety  of  problems  of  interest  by  solving  a  simple  set  of  equations. 

Briefly,  the  numerical  procedure  for  solving  an  initial-value  problem 
is  the  following.  Substituting  Eq .  (14d)  into  Eq .  (14c)  yields  a  cubic 
(algebraic)  equation  for  £ik  -  iF  vhich  may  be  solved,  at  each  z,  in  terms 
of  r^(2),  o(z)  and  rj^(z),  thus  enabling  Eqs.  (14a)  and  (14b)  to  be  stepped 
forward  in  z.  Since  in  the  absence  of  source  terms  an  input  radiation 
signal  diffracts  avay  on  the  scale  length  defined  by  the  Rayleigh  range  Zj^, 

w  r  (z)  I 

_  S  /  ^  C  \ 
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It  is  informative  to  present  the  numerical  results  with  the  distance  along 
the  viggler  measured  in  units  of  the  Rayleigh  range.  In  all  the  numerical 
results  to  oe  presented,  the  radiation  field  is  assumed  to  be  in  the  form 
of  plane  waves  at  the  entrance  to  the  viggler,  l.e.,  a(z-O)  -  0. 

Case  I 

To  begin  with.  Pig.  1  shows  the  results  for  the  following  parameters: 
beam  current,  I.  -  270  A,  r.  -  0.01  cm,  y  •  2000,  2it/k  -  10  cm,  a  -6.15 

D  DO  V  w 

1/2 

and  r  (z-0)  -  0.02  cm.  Noting  the  factor  of  2  difference  between  the 
s 

definition  of  in  Eq .  (13)  and  that  in  Ref.  4,  it  is  clear  from  Fig.  1(a) 
that  after  a  transient  oscillation  over  a  distance  of  about  20  Rayleigh 
ranges,  the  radiation  spot  size  approaches  a  value  quite  close  to  that 
obtained  with  the  two-dimensional  PEL  code  FRED  at  the  Lawrence  Livermore 
National  Laboratory  (LLNL),^  We  also  find  that  for  all  the  numerical  cases 
examined,  a  unique,  asymptotic  spot  size  is  obtained  irrespective  of  the 
initial  optical  waist.  Figure  1(b)  shows  the  spatial  evolution  of  a, 
indicating  that  it,  too,  approaches  a  constant  value  after  an  initial 
transient  behavior. 

The  solid  curve  in  Fig.  2  shows  the  evolution  of  1/e  width  of  the 
radiation  amplitude  with  a  five  mode  (m=0, 1 , 2 , 3, 4 )  source  dependent 
expansion  calculation  using  the  same  set  of  FEL  parameters.  The  radiation 
field  is  represented  by  Eq.  (3)  and  the  source  term  is  given  by  Eq .  (11). 
With  the  assumption  that  the  fundamental  mode  dominates,  only  the  6k  and  F 
of  a^  (r,z)  are  involved  in  the  source  function  and  they  are  obtained  from 
Eqs.  (14c)  and  (14d).  It  is  found  that  the  fundamental  mode  remains 
dominant  over  many  Rayleigh  lengths.  For  comparison  the  dashed  curve  in 
Fig.  2  shows  the  f u.idaoiental  mode  spot  size  of  Fig.  1(a),  and  the 
asymptotic  results  are  seen  to  differ  by  about  lOX.  This  suggests  that  the 
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single-mode  system  of  Eqs.  (1^)  may  be  regarded  as  a  reasonable  and 
accurate  simplification  of  Eq.  (3).  Henceforth  the  results  presented 
pertain  to  Eqs.  (14). 

Case  II 

Figure  3  presents  the  results  for  a  case  vhere  the  electron  beam  is 
not  matched;  l.e.,  the  envelope  of  the  electron  beam  is  modulated: 

“  “^bo  *  (16) 

vhere  ir,  is  the  amplitude  of  the  modulation  and  for  simplicity  k-  is 
b  p 

chosen  to  be  equal  to  the  betatron  vave  number^^  k^a^/ ( J 2 ) ,  neglecting 
self-fields.^^  is  the  beam  speed  along  the  viggler  axis  normalized 

to  c.  The  parameters,  typical  of  the  Advanced  Test  Accelerator  experiment 
at  LLNL,  are  I,  «  2  kA,  r,  «  0.3  cm,  y  »  100,  2jt/k  *  8  cm,  a  «  1.72, 

D  _  DO  V  V 

r^(z=0)  =  0.35  cm.  In  Fig.  J,  vhere  ^  observed  that  the 

optical  spot  size  follovs  the  modulations  in  the  electron  envelope 

apparently  identically.  Specifically,  a  number  of  cases  vere  examined  vith 

Jr, /r,  up  to  0.4.  In  all  cases  the  electron  and  optical  beams  oscillate 

vith  almost  identical  vavelength,  although  the  radiation  beam  appears  to 

lag  behind  in  phase.  Hovever,  defining  the  modulation  depth 

A  =  ((r)  -  (r)  .  ]/I(r)  +  (r)  .  J,  it  is  found  from  Fig.  3(a)  that 

max  min'  ‘  max  min'  ® 

=  0.087  vhereas,  from  Eq.  (16),  A^  «=  Although  the 

modulation  depth  of  the  electron  beam  differs  from  that  of  the  radiation 
beam,  it  is  found  that  A^  increases  in  proportion  to  Jrj^. 

More  generally,  alloving  for  the  defocusing  effect  of  self-fields, 
there  is  alvays  the  possibility  of  a  small  amplitude  ripple  on  the  electron 
beam  envelope  and  hence  on  the  radiation  beam  envelope.  In  Appendix  B, 
generation  of  sidebands  is  considered  in  a  simplified  model  and  found  to 
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have,  for  typical  cases,  an  insignificant  effect  on  the  linear  dispersion 
characteristics  of  the  primary  optical  vave,  as  implicitly  assumed  by 
employing  the  source  term  in  Eq.  (11)  in  the  present  case. 


IV.  Analysis  of  Results 

One  interesting  feature  of  the  numerical  results  is  that  in  all  cases 
the  radiation  spot  size  has  a  unique,  asymptotic  limit  Irrespective  of  the 
initial  value.  The  asymptotic  value  of  r^  and  of  o  is  determined  by  the 
fixed  points  of  Eqs.  (14a)  and  (14b);  i.e.,  at  the  fixed  point 


ck  T 

V  1 

1 

Fi  1 

^  X  2 

k  a 

-  cc 

k  a 

(k  r  ) 
vs 

(  V  oj 

R 

< 

o 

I-* 

=  0, 


ck 

4a  —  -  2 
u 


^k^a^  I  (  V  s 


k.  r  «  0. 


(17a) 


(17b) 


Combining  Eqs.  (17a)  and  (17b)  one  obtains 


(1 


(k  r  )' 
vs 


0, 


vhich,  upon  making  use  of  Eq .  (14d),  yields 


Ak 


k^r.  n 

V  b 


1/2 


U2f 


Ak/a, 


where 


and  f  =  filling  factor.  Substituting  the  expressions  for 

Ak  and  F  into  the  dispersion  relation  (14c),  one  obtains 
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a  -  (f/(3f  . 

,  ,  ,1/4  048^/2)^^'^  ,1/4.,  ,,,3/2 

(y/v) _ V  f  (l42f) 

•"s  ■  ,3/4,  ,  1/4 


2  k  Yf, 

V  B 


1/2  ,,  ,,,,,3/4  ’ 

(l43f/2) 


where  v  -  (c4^^rj^^/2c)  is  Budker's  paraaeter.  These  expressions  may  be 

used  to  obtain  the  asymptotic  spot  sire  for  a  given  filling  factor,  and 

1/2 

then  one  obtains  the  corresponding  electron  beam  radius  via  r^^  »  r^f 
To  avoid  complications  arising  at  the  outer  edges  of  the  op  cal  beam, 
where  the  field  amplitude  is  small,  typically  a  filling  factor  f  <  1/2  is 
appropriate.  It  is  also  possible  to  rearrange  the  expression  for  r^  to 
obtain 


*  (}  -  I  q]f  -  1  •  0, 


where 


1/3 


1 


143^/2 

V 


The  cubic  equation  for  f  may  be  solved  to  obtain  an  explicit  expression  for 
r^.  Noting  that  the  sum  and  the  product  of  the  three  roots  of  the  cubic 
equal  -1  and  q,  respectively,  it  follows  that  there  is  a  unique,  real  value 
for  the  asymptotic  spot  size  r^ . 

To  examine  stability,  it  is  convenient  to  define 


and  substitute  Eq.  (14d)  into  Eq.  (I4c)  to  obtain  the  local  dispersion 
relation: 
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+  2 


rck  '1 

V 


b> 


1-la 


(k  r  )‘ 
vs 


ck 

-  2  —  n 

u 


U(r^/r^)^ 
- - 


l42(r^/r^)‘ 


2’ 


(18) 


vhlch  may  be  solved  iteratively.  It  turns  out  that  for  the  parameters  of 


Case  I,  at  the  lowest  order,  the  right-hand  side  balances  the  quadratic 
term  on  the  left.  The  relevant  root,  with  tk,  T  >  0,  may  be  substituted 
into  Eq.  (14d)  to  obtain,  for  a  >  o. 


F,  -ck  ,  . 

1  V  1-io 

1  r  n  1 

1/2 

a-i 

k  a  ~  ^  /b 

VO  (kr)+(k 

V  S  V 

,2  2  L,.  2.J 

r^^)  ^(1+a  ) 

[u(U«2)1/2]1/2 

(k  r  )^ 

'vs 

(19) 

(k  r  )^  +  2(k  r,  )^ 
'vs  V  b 

Perturbing  Eqs.  (14a)  and  (14b)  about  the  fixed  point  and  making  use  of 
Eq.  (19),  it  is  simple  to  show  that  the  perturbation  is  spatially  damped, 
thus  indicating  the  stability  of  the  fixed  point.  The  algebraic  details 
are  relegated  to  Appendix  C. 

Another  aspect  of  the  results  which  is  of  interest  pertains  to  the 

nature  of  the  phase  fronts  and  the  flux  of  optical  power  in  the  asymptotic 

region.  From  Eqs.  (7)  and  (8)  it  is  simple  to  check  that,  in  differential 

2 

form,  the  surfaces  of  constant  phase  are  given  by  (w/c  ♦  Ak)6z  (2ra/r^)6r 

=  0,  and  hence,  noting  that  Ak,  a  >  0,  the  wavefronts  are  divergent  in  the 

direction  of  propagation.  Consistent  with  this,  there  is  a  nonvanishing 

transverse  component  of  the  Poynting  flux.  Specifically,  for  <  1  the 

ratio  of  flux  of  optical  energy  in  the  transverse  direction  to  that  along 

2 

the  2  axis  is  -  or/kr^  <<  1. 
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V. 


Conclusion 


Based  on  the  results  presented  herein,  the  simplicity  ajid  accuracy  of 
the  single-Dode  Gaussian-Laguerre  approximation  to  the  solution  of 
Kaxvell's  equations  has  been  demonstrated.  It  is  shovn  that,  in  the 
exponential  gain  regime  of  operation  of  an  PEL  amplifier,  there  is  a 
unique,  asymptotic  spot  size  for  the  radiation  beam  irrespective  of  that  at 
the  entrance  of  the  wiggler.  There  is,  however,  a  transverse  flux  of 
optical  power.  It  is  shown  analytically  that  the  asymptotic  profile 
(i.e.,  the  radius  and  the  curvature  at  large  z)  is  stable  to  small 
amplitude  perturbations.  Vith  a  spatially  modulated  electron  beam 
envelope,  that  of  the  optical  bean  is  found  to  oscillate  on  the  sane 
spatial  scale. 
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Appendix  A:  Source  Ter« 


In  this  appendix,  the  details  of  the  evaluation  of  the  source  term  S 
in  Eq.  (11)  are  presented. 

The  PEL  source  current,  J^(r,e,2),  In  a  linear  viggler  is  given  by 
Jjj(r,e,z)  -  -  |e  e,z)v^ 


.2-  -ik  z 

I  '' 

2y®  c 
o 


where  electron  self-field  effects  are  neglected.  Taking  the  convective 
time  derivative  of  Er .  (A2),  and  incorporating  the  linearized  version  of 
Eq .  (A3),  one  car.  arrive  at  the  following  equation  for  the  perturbed  beam 
densi ty , 


14 


d^in 

- 

D 

dt^ 

■ 

c 

dz 

1  C  4. 

Uz  *  2  at 

c 

J  pond 

(A4) 

vhere 

*pond 

-  I'lv 

.  2 
Arm  c 

0 

e 

i  [(k 

*  k^)z 

-  wt  1 

J  +  c.c. 

With 

the  assumption 

that 

A(r,e, 

z)  is  a  slovly  varying 

function  of 

i.e. ,  j31nA/3z |  <<  k^ 

« 

k,  Eq.  (A4) 

becomes 

d^&Hb 

i 

[fk  .  k, 

D  V 

k 

V 

,ke 

L  V 

^  J  J 

(A5) 

dt^ 

'-222 
2y  m  c 
o 

*  C.C* 

vhere  the  resonance  condition,  u  v^(k.  +  k^)  is  used. 

For  a  near  Gaussian  radiation  field  in  the  exponential  gain  regime, 


A(r, e,z)=A^(r,e,z)=A^(o)exp< 


(  Ak(Zj)-ir(z^) ]dz^- 


1 l-io(z) J 


and  assuming  dk,  F,  «  and  r^  are  slovly  varying  functions  of  z,  Eq .  (A5) 
can  be  integrated  immediately  to  give 


inb 


jel^n.  A  A  k  k  i  |(k  +  k  )z  -w 

'  '  b  V _ V  L  v' 

2  2  4  2  ® 

2y  m^j  c  (Ak-iF) 


c .  c . 


(A6) 


When  Eq.  (A6)  is  substituted  into  Eq.  (Al),  taking  into  account  the 
usual  difference  of  Bessel  functions  for  a  planar  viggler,  and  Eq .  (10)  for 
the  beam  profile,  the  source  function  in  Eq .  (1)  is  then  given  by  Eq .  (11). 


15 


Appendix  B:  Sideband  Generation 

In  this  appendix  generation  of  sidebands  to  the  primary  optical  vave, 

due  to  the  spatial  modulation  of  the  electron  beam,  is  analyzed.  It  is  to 

be  emphasized  that  the  folloving  analysis  is  intended  merely  to  show  that 

the  dispersion  characteristics  of  the  primary  optical  wave  are  only 

2 

slightly  modified  (-  (SN^/N^)  ]  for  typical  experimental  parameters,  as 

implicitly  assumed  in  applying  the  results  of  Appendix  A  to  the  case  of  a 

modulated  electron  beam  in  Section  III. 

The  development  of  the  linear  theory  herein  generalizes  that  of 
12 

Sprangle  et  al.,  to  which  reference  should  be  made  for  further  details. 

The  form  of  the  vector  potential  of  a  planar  viggler  employed  in  this 
appendix  is  slightly  different  to  that  given  by  Eq .  (12): 

X’ 

where  is  purely  imaginary,  and  that  of  the  linearly  polarized  radiation 
field  is  taken  to  be  of  the  form 


A  =  |a  exp  ik  z-ico  tj  +  A  exp^ik  z-imt 


t  exp  ik  z-i(i)  t  *  c.  c.  >  e  , 
O  \  /  I  ^ 


where  it  is  assumed  that  the  electron  density,  modulated  at  the  betatron 

wavelength  2n/k:-,  has  the  simple  form 
p 

iN 


“  r  1 

n^  =  exp(ikgZ)  +  c.c.J, 


with  k^  <<  k  <<  k,  and  k  »  k  -►  k.,  k  =  k  -  k.. 

gw  ♦  6  -  g 
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Folloving  Ref.  12,  tl»e  vave  equation  is  found  to  be 


cr  :  3^  .  4n|er  ,  , 

Irr  -  "2  rr  ■  — 2J  -  ‘  ”2 

iz  c  3t  Y  c 

o  00 

vhere  r  is  the  relativistic  factor  in  the  absence  of  the  radiation  field, 
0 

2  1/2 

-  (4 Jin^  |e  I  /b^)  ,  and  in  is  the  density  perturbation  caused  by  the 

radiation.  Note  that  the  velocity  v^^  along  the  vlggler  axis  is  not 

affected  by  the  betatron  oscillation  and  hence  y  .  to  lowest  order  in 

o 

2  2 

leA^/Y^B^c^ I  ,  is  not  a  function  of  z.  Defining  the  ponderomot i ve 
potential 


Y  m  c 
o  0 


1 


the  momentum,  continuity  and  Poisson's  equations  may  be  combined  to  obtain 


.2  V  .3n  .  ,  4itn  le  •  11  ,3n  .  .. 

-  - - ^  ‘  ^  5r 

Qt  0 

o  O  Z  0  o  z 


-jel  l_  !zo  3_'i. 

m  Y  3z  ^0  Uz  *  2  3tj  pond' 


2  2  -1/2 

where  y^  =  scalar  potential,  and  terms  such  as 

2  2  2 
3  n  /3z  ,  which  are  on  the  order  of  k_,  have  been  neglected, 
o  p 

Writing  .  k  *  kg,  k_  =  k  -  kg, 
in  =  exp  l^i  ^k^  +  k^jz-iu  tj  +  in_  exp  |^i  ^k  +  k^jz-iw  tj 


in^exp  [i  (k  . 

^w]^  i“t] 


noting  that,  on  the  left-hand  side  of  Eq .  (Bl),  the  ratio  of  the  fourth  to 

the  third  term  is  on  the  order  of  k_/k  <<  1,  one  finds  that 

P 
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2 

vhere  e  «  (5N^/2N^),  and  then  simple  to  shov  that, 

2 

correct  to  0(c  ),  the  dispersion  relation  is  given  by 


vhere 


is  the  usual  matrix  element  for  the  primary  vave,  m, ,  =  m„_(k  ), 

11  22  + 

“33  "  “22^^-^' 

“12  "  ®12^^+’ 
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to  -  [k  +  k  Iv 

“bo 

g  ro 

t  v)  zo 

♦  2  f 

Y  Y  / 
o  z 

“bo 


2  1/2 

and  to,  -  N  /m  )  .  Note  that  vith  the  definition  chosen  for  A 

^o  0  o 

this  appendix,  a^  -  (eA^/m^c^)^  <  0. 

To  proceed  along  the  lines  of  Ref.  12,  it  is  convenient  to  vrite 

®22  '  ^22  *  ^22’ 


vhere 


’22 


Ico  -  |k  +  k.  Iv 

v)  zoj 


2  2 
“bo 


,2  0)^  “bo 

2  (  {^  -  ~2  *  2 
r^y^)  c  Y^c 


). 


and 


'22 


is  the  "coupling"  term.  The  dispersion  relation  then  becomes 


^22  dispersion  relation  for  uncoupled  electromagnetic  and  space 

charge  waves.  The  right-hand  side  of  Eq.  (B2)  Introduces  the  PEL 
Interaction  and  coupling  to  sidebands,  and  Its  effect  is  Included 
iteratively.  At  the  lowest  order,  -  0  for  some  (w.k).  Substituting  In 
the  right-hand  side,  the  second  set  of  terms  vanishes;  the  term 
proportional  to  C22  survives. 

Substantial  modification  of  this  dispersion  relation  is  expected  if 


1 


/ 

0 

2N 

oj 

^33 


<<  1 


1 . e . ,  if 


-  2 


-3/2 

o 


1/2  3/4 

Y,  Y^ 


For  typical  experimental  parameters,  the  right-hand  side  of  this  equation 
exceeds  unity,  whereas  «  1,  implying  the  insignificance  of  the 

effect  of  modulation  on  the  dispersion  relation. 
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Appendix  C:  Stability  Analysis 


The  purpose  of  this  appendix  is  to  establish  the  stability  of  the 
fixed  point  (r^,«)  of  Eqs.  (14). 

Perturbing  Eqs.  (14a)  and  (14b)  about  the  fixed  point  and  making  use 
of  Eq .  (19),  it  is  seen  that  the  perturbation  evolves  according  to: 


So) 


d(k^z)  KLk 


2  ^11  ^12 
®21  ®22 


)  (a. 


where  x  -  (k^r^)^,  y  = 


-  a 


(ckv/“) 


‘11 


12 


x-y 


3o 


oJCil, 


-(Uo^)  ^ck^/'i)]y(2x+y) 


2  2 

X  (x+y) 


3x 


_  _  i_  Y 

®21  =  o  3a 
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a[ck^/u)] 

(x-^y)^ 


fc(-i) 


and 


1/2 


•  fl  (^ 

0-1  1+  11+0  J 


2 (!.«')] 

Assuming  that  So,  6x  -  exp(Xk  z),  one  finds  that 


P  l/2.,l/2  x*2y  tx*y 

[l.(u«2]  ] 


,  .1/2 


X  =  - 


a  ^ck^/ojj  (x  +  2y ) 


(x*y)‘ 


♦52*5^ 


(x-y) 


2  ^  ^2"  ^1 


,2  1/2 

-  sA  , 
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vhere 


S 


1 


ax 

®  3^ 


H 

ck  /(o 

1  y(y*2x)  g 

. 

2/  x2 

X  (x+y) 

*  ax 

(Note  that  all  the  variables  in  this  appendix  are  evaluated  at  the  fixed 
point.)  Making  use  of  Eq.  (Cl)  it  is  simple  to  show  that  X^/a  -  ax^/aa  >  0, 
a(Xj^  -  oiXj)/ax  <  0,  whence  >  0  and  hence,  noting  that  >  0,  and 

that  the  perturbation  solution  for  Eq .  (18)  implies  <  a(ck^/a))/(x+y ) ,  one 
finds  that  Re  X  <  0,  thus  indicating  the  stability  of  the  fixed  point  to 
small  amplitude  perturbations. 
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PigTjre  Captions 


Fig.  1 


Fig.  2 


Fig.  3 


Spot  size  (fj)!  phase  shift  (&k),  and  gain  (F)  vs.  distance 

along  the  viggler.  z  is  normalized  to  the  Rayleigh  range  z  .  In 

K 

(c)  and  (d)  the  number  on  the  ordinate  must  be  multiplied  by  10~^ 
and  10  ^  to  obtain  the  actual  value  for  &k/k  and  F/k.  , 

V  V 

respectively. 

(l/e)-vidth  of  the  optical  field  vs.  distance  along  the  viggler. 
Solid  curve:  5  mode  system;  dashed  curve:  1  mode  system. 

Spot  size  (tg))  <*1  phase  shift  (Ak),  gain  (F),  and  radius  of 
electron  beam  (r^^)  vs.  distance  along  viggler. 
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